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Abstract 



This paper analyzes the scattering theory for periodic tight-binding Hamiltonians perturbed 
by a finite range impurity. The classical energy gradient flow is used to construct a conjugate 
CN \ (or dilation) operator to the unperturbed Hamiltonian. For dimension d > 3 the wave operator 

is given by an explicit formula in terms of this dilation operator, the free resolvent and the 
■ perturbation. From this formula the scattering and time delay operators can be read off. Using 

the index theorem approach, a Levinson theorem is proved which also holds in presence of 
embedded eigenvalues and threshold singularities. 

1 Introduction 

> \ The purpose of this work is to present the scattering theory for a quantum particle described by a 
["J- " tight-binding Hamiltonian H = Hq + V acting on the Hilbert space £ 2 (Z d ) where H is a periodic 
tJ - ■ operator with a single band and V is a finite rank perturbation. Most of the present work is focusing 
on the case of dimension d > 3. All along this paper, it will be assumed that the Fourier transform of 
Hq acts on L 2 (T d ) as a multiplication operator by a real analytic Morse function £(k) having only one 
maximum and one minimum. Operators of this type appear in solid state physics as effective one- 
band Hamiltonians for electrons or holes in periodic media. The analyticity reflects the exponential 
decay of the hopping amplitude of the particle and the Morse condition is generic. 

X: 

1.1 Main Results 

Scattering theory for a Schrodinger operator with a periodic potential has already been considered 
|New2| IBY| IGNj . The above scattering problem has also been addressed in the physics literature, 
for example in |Eco] . The present work is going further. Initially, it was motivated by the remark by 
Kellendonk and Richard [KR1] that Levinson's theorem |Lev] relating the number of bound states 
to the total scattering phase can be interpreted as a special case of the Atyiah-Singer index theorem. 
As it turns out, this nice basic idea requires a substantial amount of technicalities when it comes 
to mathematical justification |KR2| IKR3[ IKR4] . For indeed, the global character of Levinson's 
theorem requires several technical steps. First, a complete control on the nature of the singularities 
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of the Green function of Ho is needed, a task that is easy on the continuum, but more involved in 
the present case. In addition, it requires a conjugate operator to H in order to shift the energy, 
replacing the dilation operator used for the continuum situation. Moreover, the potential term V may 
create embedded eigenvalues and threshold singularities that must be analyzed thoroughly since they 
contribute to Levinson's theorem. At last, for the index theorem to apply, it is necessary to prove that 
both the wave operator and the scattering matrix can be expressed as suitable continuous function 
of the energy and dilation operators. This was achieved in [KR3J through an explicit calculation in 
dimension d — 1 using the techniques of [Jen] . In higher dimension |KR4] . the explicit calculation 
turns out to be harder, but it is possible to prove sufficient regularity of the wave operators and the 
scattering matrix. Unfortunately, the work [BY] is insufficient to implement this program completely 
for the class of models considered here. The present paper is supplementing these points. 

In the light of the previous introduction, the main results of the paper can be summarized as 
follows: 

• A lattice analog of the dilation operator is constructed (see Theorem [T]). It is a self-adjoint 
unbounded operator A such that i[Hq, A] = F(Hq) where F is a positive function on the 
spectrum of H vanishing only at the band edges. 

• Explicit expressions for the wave operators, the scattering matrix and the time delay operators 
are obtained in terms of Ho, V and A (see Theorem [21 Theorem [3] and Theorem H]). 

• A series of results concerning the existence of embedded eigenvalues and threshold singularities. 

• A Levinson type theorem is derived which is now briefly described (see Theorem [5] for details). 
The essential spectrum of H is the spectrum [E_,E + ] of Ho and E± are called the band 
edges or thresholds. Let P pp be the eigenprojection on the eigenvalues N be the total number 
of eigenvalues, including the embedded and the threshold eigenvalues (the latter have to be 
distinguished from threshold resonances). Further let S be the scattering operator and let 
T = — iS^ 1 ^, S] be the time delay operator seen as acting on £ 2 (Z d ). Finally let m± e {0, 1} 
be the degeneracies of the threshold resonances at E± also called half-bound states singularities 
(higher degeneracies are possible, but not dealt with here). Then for d = 3 and isotropic 
extrema of S, 

±Tr{S- l [A,S\) + Tr(P pp ) = - ^ - ^ . (1) 

For d > 5 it is proved that always m± = and that (CEJ) holds. For d — 4, (CD) is proved under 
the hypothesis that m± = 0. 

As already pointed out, this paper is restricted to dimension d > 3. Dimensions d — 1 and d = 2 
require a detailed asymptotic expansions of the free Green function near the band edges [Newl[ 
IBGDWl IKR3] . The one-dimensional case has been treated in |CK[ IHKSj . The case d = 2 will be 
addressed in a future publication. 
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While most of these results are technically new, similar results have been already obtained in 
the past. The scattering problem in M. d with standard Laplacian perturbed by a decaying potential 
together with a proof of Levinson's theorem can be found in standard textbooks such as |Newll |R5) . 
In this situation, also threshold resonances have been analyzed in details (see [Bol j for a review). The 
scattering of an electron in a periodic potential by an impurity has been considered by physicists for 
a long time, in connection with the transport properties of semiconductiors. This theory is based 
on the KKR equations (for Korriga, Kohn and Rostoker) and in this context Levinson's theorem 
given in equation (TjQ) above is also known under the name of Friedel sum rule. This was investigated 
by Newton |New2] for d = 3 with an impurity potential that could lead to threshold resonances, 
but not embedded eigenvalues. Levinson's theorem for periodic potentials in dimension d — 1 was 
proved by Firsova |Fir] . The mathematical aspects of scattering theory in a periodic potential has 
been considered by Birman and Yafaev [BY] and bears many similarities with the present approach. 
However, the latter work does not deal with the critical points of the band functions, and hence 
does not lead to a proof of Levinson's theorem. The scattering of a lattice electron by a localized 
impurity is also addressed in the book by Economou [Ecoj . however, the explicit formulas obtained 
in the present work for the scattering matrix and the wave operators are missing. The problem 
of embedded eigenvalues has long been considered as an irrelevant curiosity. Indeed, even though 
they are non-generic within the class of finite rank perturbations considered here (as follows from 
the arguments in Section r3.2l) . they may occur in practical devices, in particular, when a compact 
part of the lattice is inaccessible to a particle coming from the outside (see Example H] below). The 
eigenvalues of the Hamiltonian inside have an influence on the scattering outside in various ways, as 
can be seen in equation (CTJ. An analogous effect occurs for microwaves reflected by a cavity, as was 
shown, for instance, in [DSFJ. 

1.2 Strategy of proofs 

As suggested by the formula i[Hq, A] = F(H ), the conjugate operator A is the generator of an energy 
shift. For its construction the classical energy gradient flow is slowed down near the band edges which 
are also called thresholds. This flow can be implemented as a strongly continuous one-parameter 
group of unitary operators in the Hilbert space and A is then simply the generator of this group. The 
unitary implementation of a vector field has already been carried out in [HS,ABG], however, these 
constructions excluded energy surfaces with critical points. Removing this constraint is crucial for 
the proof of Levinson's theorem and this is probably the main conceptual contribution of this paper 
to the scattering theory in periodic lattices. The proof also covers dimension d = 2. 

The introduction of the conjugate operator A is closely linked to an important tool of calculation 
used in this paper, namely the reseated energy and Fermi surface (REF) representation giving an ad- 
equate spectral representation of both Hq and A. It shows that the Hilbert space £ 2 (Z d ) is isomorphic 
to L 2 (1R) <g> L 2 (S, v) where R is a rescaled energy variable and E is a reference Fermi surface given 
by some level set of £ furnished with a Riemannian volume v. The influence of the critical points on 
the dynamics defined by A lies on a set of zero Lebesgue measure explaining why the unitary group 
e ltA is globally defined. While Morse's theory shows that the topology of the level sets changes at 
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the passage through a critical value, the previous result, on the opposite, shows that the topology 
of the Fermi surface does not play any role for the Hilbert space isomorphism. Changing Hq into 
B = f(H ), for a suitable function / such that i[B, A] = 1, leads to a representation where B is the 
multiplication by a variable 6 6 M, called the rescaled energy while A becomes the derivative — idf, 
and acts as an infinitesimal rescaled energy shift. 

Given the REF representation, it is possible to compute all standard objects of scattering theory 
explicitly. In order to limit the technical difficulties, this work will be restricted to dimension d > 3 
and to a compactly supported perturbation. The wave operator, minus the identity, is then an 
explicit continuous function in A and B with values in the algebra of compact operators on L 2 (£, v). 
From this formula and the invariance principle, an expression for the on-shell scattering matrix is 
then readily deduced. Up to an explicit partial isometry, it is a finite dimensional unitary matrix 
expressed in terms of the perturbation V and the Green function of Hq. The spectral property of 
the time delay operator linking it to the resolvent also follows form this analysis. This allows to 
give a first short proof of Levinson's theorem by a contour integration argument when there are no 
embedded eigenvalues and no threshold singularities. 

However, both threshold singularities and embedded eigenvalues may occur for adequate choices 
of V, sometimes with physical meaning. This situation is covered by the second proof of Levinson's 
theorem which follows closely the .fT-theoretic arguments of Kellendonk and Richard [KRlj IKR2j . 
Following these authors, a C*-algebra £ is generated by continuous functions of A and of B with 
values in the compact operators on L 2 (S, v) and having well-defined limits at ±oo which coincide in 
the four corners A = ±oo and B = ±oo. It contains the ideal J of those functions vanishing at oo 
and the extension is precisely by the C*-algebra A of operators fibered over A or B, again coinciding 
in the four corners. A large amount of effort is then dedicated to proving that the wave operator 
belongs to the Toeplitz extension £, even when embedded eigenvalues and threshold singularities are 
present. It follows that the wave operator is a lift of the scattering operator which combined with 
contributions stemming from the thresholds is an element of A. This leads to a .ff -theoretic version 
of the proof of Levinson's theorem. 

Notations: As usual, \A\ = (A*A)^, $teA= \{A + A*) and ^smA = ±(A-A*) for any operator A. 
Furthermore, throughout there is a rescaled energy variable b = f(E) associated with the bijection 
/ from the spectrum of the unperturbed operator to R which is defined in (jSj) below. For objects 
depending on energy both E and b will be used as indices, for example = Pe, ILj = He, Cf, = Ce 
and so on. 

Acknowledgments: We thank A. Knauf, S. Golenia, S. Richard and J. Kellendonk for numerous 
comments. The work of J. B. was support in part by NSF Grant No. 0600956 and 0901514, that 
of H. S.-B. in part by the DFG. After the first version of this paper was submitted, two papers on 
related matters appeared on the archives. The article [BPLj provides a variational principle giving 
an upper bound on the number of eigenvalues (including embedded ones) and [KKNj analyzes the 
Friedel sum rule via the spectral shift function. 
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2 Analysis of the unperturbed lattice Hamiltonian 



2.1 Unperturbed Hamiltonian and its energy band 

The tight-binding Hamiltonians considered in this work act on the Hilbert space £ 2 (Z d ) of square 
summable sequences of complex numbers indexed by the <i-dimensional lattice Z d . The free Hamil- 
tonian Hq on £ 2 (Z d ) is supposed to be of the form 



where the £ n 's are the Fourier coefficients of a real-analytic real- valued function £{k) = J2nez d elkn £n 
on the o?-dimensional torus T d = IR d /(27rZ d ). Hence we restrict ourselves to a free operator with 
a single band. As H is translation invariant, it is diagonalized by the discrete Fourier transform 
T : £ 2 (Z d ) — > L 2 (T d ), where L 2 (T d ) is the Hilbert space of square integrable functions on T d . It is 
densely defined by 



and is unitary. Now Hq = FHqF* is a multiplication operator on L 2 (T d ) by the function £. The 
main hypothesis on Hq are expressed in terms of this function £. The set of critical points S* C T d 
at which the gradient V£ w.r.t. the euclidean metric vanishes is finite due to the analyticity of £ 
and each critical point is supposed to be non-degenerate, namely for any k* e S* the Hessian £"{k*) 
is a real symmetric invertible al x al matrix. In other words, £ is supposed to be a so-called Morse 
function |Nic] . Recall that the index of a critical point k* is the number of negative eigenvalues of 
£"(k*). Then the Morse inequalities state that the number of critical points with index p is larger 
than or equal to the Betti number /3 P of the torus T d , which is equal to the binomial coefficient al over 
p. In particular, there must exist critical points of £ with signature p for every p — 0, . . . , d. For the 
discrete Laplacian, the energy band £{k) = 2 J^ =1 cos(kj) is a Morse function for which the Morse 
inequalities become equalities. We also assume that there are only two critical points k*_ and k* + of 
definite signature corresponding to the minimal and maximal values £L = £{k*_) and E + = £{k* + ) 
of £. Hence all other critical points k* e S* are supposed to have critical values £(k*) in (E_,E + ) 
and to be of indefinite signature. Note that £(S*) is the set of all critical values. A spectral interval 
is called non-critical if it does not contain any critical value. 

2.2 The classical energy flow 

Let F : [E_,E + ] — > IR> be a real analytic function vanishing only at the band edges E_ and E + 
and satisfying F(E - EJ) < C\E - E_\ and F(E + - E) < C\E + - E\ for some constant C. Below 
we will choose 





(2) 



m&L d 




but this particular choice will only become relevant for the calculation of the wave operators in 
Section 13.41 Then let X be the vector field on T d defined by 

m = F(S(k)) ™j£L , ke T d . (4) 

Apart from the factor F o £, the vector field X is precisely the one used in the standard argument 
of Morse theory |Nic] as well as in the proof of the coarea formula [SakJ. As £ and F are smooth, 
this vector field is smooth away from the set §* of critical points. At the critical points k± with 
extremal energy £{k±) = E±, the function k y F(£(k± + k)) vanishes linearly by the assumption 
on F and hence the vector field has a source or a sink there. At all other critical points with critical 
values lying inside the band [E_, E + ], the vector field X has a singularity which has to be dealt with 
below. Let 9 b : T d \ S* -> T d be the flow of X, that is, d b 9 b = X o 6 b and 9 = id. The somewhat 
unconventional choice of b as notation for the time parameter is due to its interpretation as rescaled 
energy variable below, which is dual to the spectral parameter a of the dilation operator A. The flow 
9 b is not complete because an orbit can reach one of the critical points with indefinite signature in a 
finite time. Choosing orbits which stay away from these critical points or times which are sufficiently 
small, one can calculate the flow of energy along the orbits. By the definition of the vector field X, 

d b £{9 b (k)) = F(£(9 b (k))) . 

This equation shows that the flow 9 b maps constant energy surfaces to constant energy surfaces. 
Moreover, the energy flow is governed by a simple ordinary differential equation of first order which 
can be integrated. Choosing some reference energy E r G E + ), it leads to the following invertible 
function 

f(E) = / T77~T • (5) 



F(e) 



Then b = f{£{9 b {k))) - f(£(k)) and 



£(9 b (k)) = f- l (b + f(£(k))). (6) 
If F is given by equation (j3J) and if E r = (E + + EJ)/2, it gives 

^HK)' /-<»)=* + * *»"<»). F ^(6)) = -^, (7) 

where A = (E + — EJ)/2. By restricting 9 b to an adequate subset of T d , a complete flow can be 
constructed. Let S be the union of S* and of the set of points reaching one of the critical points 
k* G S* in finite time (either positive or negative). It is important to remark that, under this flow, 
almost all points reach the maximum and the minimum eventually, but it takes an infinite time to 
do so. Therefore the finite time condition is a strong constraint. In fact, S is the union of S* and the 
stable and unstable manifolds of all critical points of indefinite signature. 
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Proposition 1 The set S is compact and has zero Lebesgue measure. The flow 9j, : T d \ S — > T d \ S 
is defined for all b 6 M., that is, X is complete on T d \ S. In addition, lim^ioo $b(k) = k± for all 
k G T d \ S. Furthermore, for any open neighborhood U of S there exists an open subset V C U which 
contains S \ {A;* , k* + } and is invariant under the flow 9. 

Sketch of a proof. The vector field X is gradient-like in the terminology of [Nic] (it is actually a 
gradient vector field). Hence |Nicl Section 2.4] shows that lim^-too Ob{k) e S* and that the stable and 
unstable manifolds of all critical points of indefinite signature are locally smooth submanifolds of T d . 
For each critical point, the sum of the dimensions of the stable and unstable manifolds is equal to d. 
Along the flow on these submanifolds the energy increases with a finite speed, except in neighborhoods 
of k±. Hence either the submanifolds reach another critical point in a finite time (non-generic) or the 
points k* ± in infinite time. Consequently, the points k* ± compactify the stable and unstable manifolds. 
As the number of critical points is finite, the set S is compact with zero Lebesgue measure. To prove 
the last statement of the proposition, let k* be a critical point of indefinite signature. Then let V(k*) 
be an open neighborhood of k* contained in U. Then V = {J k * U&gr ^b(V(k*)) is an open set that is 
invariant by the flow. A compactness argument can be used to show that V C U by choosing V(k*) 
sufficiently small. □ 

The level set of £ corresponding to an energy E e (E_, E + ) is defined by 



S E = jfcGT^S £{k) = E} 



These level sets will be called the quasi-Fermi surfaces. This terminology is introduced to stress 
that He is a strict subset of Fermi surface £~ 1 (E) because the points on the stable and unstable 
manifolds of all critical points with indefinite signature are excluded. However, the difference is 
only of measure zero. A reference quasi-Fermi surface will be taken at energy E r and denoted by 
E = T, Er . Because the singularities are excluded, the sets T, E are smooth open submanifolds of 
T d of codimension 1 which, for d > 2, have several connected components. Now the flow 6^ maps 
these connected components diffeomorphically into each other. By the above arguments, for each 
energy E, there is a time b = f{E) such that the flow 6b maps the reference quasi-Fermi surface S 
diffeomorphically into Eg. Consequently we have: 

Proposition 2 For E e (E_,E + ), the map Of(E) '■ E — > E^ is a diffeomorphism. 

For our purposes below, we will also need properties of the divergence of A. A straightforward 
calculation gives 

div(A")(fc) = F'(S(k)) + F(S(k)) (AI£L_ 2 V£(kW { WE(k)) 



\V£(k)\ 2 |V£(A;)| 4 

where Dirac notation is also used for vectors in M. d . Near a critical point k*, one has V£(k* + k) 
£"{k*) k + 0{k 2 ) } leading to 

div(X)(*. + *) = F'(£(k* + *)) + F m - + *)) - 2 + ( W -')) 
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Both F and F' are regular, hence \div(X)(k* + k)\ < C/\k\ 2 and thus div(X) is an integrable function 
for dimension d > 3. Furthermore, near the extrema k±, namely at the band edges, £{k± + k) = 
£(k* ± ) + ^(k\£"(k* ± )\k) + 0(\k\ 3 ), F'(£(k* ± + k)) = T2 + 0(\k\ 2 ) and F(£(k* ± + k)) = T{k\£"(k* ± )\k) + 
0(\k\ 3 ). Therefore, setting 

Tr{£"{k*±)) (k\£"(k* ± )\k) (k\£"(k* ± ) 2 \k) - 2 (k\£"(k±)\k) (k\£"{k*±) 3 \k) 

9± ^ ' ~ (k\£"(k*±y\k) 2 ' 

leads to 

div(X){k* ± + k) = T (2 + g ± (k) + 0(\k\)) . (8) 

The functions g± are homogeneous of degree and can thus be seen as functions on the sphere §> d ~ l . 
In dimension d — 1, one has g±(k) = —1. In higher dimension, g±(k) = d — 2 whenever £' f (k±) is a 
multiple of the identity (isotropy of the extrema). Otherwise g± are non-trivial. 



2.3 Construction of the dilation operator 

The aim of this section is the construction of an unbounded conjugate (or dilation) operator A such 
that i[A, Hq] = F(H ) where F is as above. The basic idea is to implement the flow 6^ of X in 
L 2 (T d ) as a strongly continuous group of unitaries. Let T> denote the set of smooth functions on T d 
vanishing in some neighborhood of S. Since S has zero Lebesgue measure and is compact, T> is dense 
in L 2 {T d ). Furthermore, Proposition [1] implies that every function in T> vanishes on a flow invariant 
open subset containing S \ {k*_ , k* + }. Hence for e T>, the following operator can be defined 

{W b (j>){k) = ex pQ^ dudw(X)(9 u (k))^ <j>{O b {k)), (9) 

because the singularities of X are not reached, due to the restriction on the support of 0. The 
unitarity of Wb follows from the change of variables k i— y 6b{k) and from the Jacobian formula 

det{e' b {k)) = ex P^ du div{X){6 u {k))) . (10) 

This latter relation follows from integrating d& In det = div(X)(9b(k)) with the initial condition 
det(^o) — 1- Furthermore, the group property 0b o 6 U = 0b+ u immediately implies WbW u = Wb+ U - It 
can be checked, by a direct calculation, that ||W&0|| = \\<f>\\ for G V. In addition, using the Lebesgue 
dominated convergence theorem, lim^o VVb0 = for e T>. It follows, from a 3e argument, that 
Wb can be extended as a one-parameter, strongly continuous group of unitary operators on L 2 (T d ). 
By Stone's theorem the generator A = ^dbW b \ b =o is self-adjoint and W b = exp(ibA). Also [BR, 
Corollary 3.1.7] implies that T> is a core for A because T> is left invariant under W&. The derivation 
of equation ([9]) leads to 



8 



A<J> = - \x(<j>) + - div(X) <j>j , (11) 

where X((j)) = (X|V)0 is the action of the vector field on the function e V. Note that the 
multiplicative (zero order) operator | div(X) is needed to make the r.h.s. of (fTTT) symmetric w.r.t. 
the scalar product in L 2 (T d ). The desired commutator property i[A, Ho] = F(Hq) now follows 
directly from ( TTTl) because i[A, H Q ] = X(£) = F(H ). This can be summarized as follows: 

Theorem 1 Let £ be a Morse function with only one maximum and one local minimum and let F 
be a smooth function vanishing linearly at the two extremal values £L and E + and nowhere else. Let 
Wb be defined by (Q for <fi e T> and with X and Of, given by (J4]) and its flow. Then Wb is a strongly 
continuous one-parameter group of unitary operators on L 2 (T d ). Its generator A = ^dbWb\b=o is 
self-adjoint with core T> and satisfies 

i[A,H ] = F(H ), i[AJ(H Q )} = 1. 

A few comments conclude this section. The vector field X defined by (jlj) has singularities stem- 
ming from critical points of £ with indefinite signature (where F does not vanish). This leads to 
singularities in both the principal and subprincipal symbol of the differential operator A as given in 
(TTTj) , As shown in Section |2^2| the singularity of the principal symbol is integrable in dimension d > 2 
while the subprincipal symbol is integrable for d > 3. It has been shown above that this does not 
prevent (TTTT) from defining a self-adjoint operator. There is another similarity between A = T*AT 
and the usual dilation operator used for the Laplacian in L 2 (IR d ). Let Xj = J r *XjJ r be the operator 
on £ 2 (Z d ) associated with the jth component Xj of X. Also let Q = (Qi, . . . ,Qd) be the position 
operator defined by Qj 4>(n) = rij <f>(ri), for n e 7L d and decreasing sufficiently fast. Then the Fourier 
transform of the r.h.s. of ( TTTT) leads to 

1 d 

^ = ^(XjQj + QiXi). (12) 

3=1 

Comparing with the usual dilation operator on M. d , Xj can be interpreted as the lattice analog of the 
jth component of the momentum operator. 



2.4 Change of variables and REF representation 

This section is devoted to the definition and the properties of the rescaled energy and Fermi surface 
(REF) representation. The proof of Theorem [1] was mainly based on the change of variables 6b ■ 
T d T d with Jacobian (fTQ|) . It will be supplemented by the coarea formula (see e.g. |Sak] for a 
proof and note that S is of zero measure). If ve denotes the Riemannian volume measure on 
(induced by the euclidean metric on T d ), 

[ dk<P(k) = [ E+ dE [ v E {da) —1— <f>(a) . (13) 
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This holds for in the set D. For the reference energy surface £ = E^ r , the measure is simply 
denoted by v = v Er . The coarea formula leads to the following: 



Lemma 1 Let <p G V. Then its integral can be written in the following three equivalent ways: 



dk (f)(k) 



db / v{da) 



det(0i 



X{6 b {a)) 4>{6 h {a)} 



db / v(da) exp ( / du div(X)(6 u (cr)) 



X(a) 



<f>(0 b (o-)) 



E- 



dE J v{da) 



det(0'. 



HE) 



T CT S 



\VS(9 f(E) (a))\ 



<t>{e m {°)) , 



(14) 
(15) 
(16) 



where ^|r CT s denotes the derivative of 9 b restricted to the tangent space of '£ at a (so that this is a 
(d — 1) x (d — 1) matrix). 



Proof: Starting from the coarea formula (fT3l) . the substitution b = f(E) given in ([5]) and the 
diffeomorphism of Proposition [2] will be used in the following change of variables: 



dk 4>{k) 



jd 



db 
db 



F(f-\b)) 

\V£(6 b (a))\ 



f-Ht>) 
u(da) det(^| Ta s 



{a)). 



In the second equality, the identity F(f 1 (b)) = F(£(a)) for o G Hf-i(p) was used. Replacing 
the definition of X already shows ( IT41 as well as ( JT6l) . Next 6' b can be decomposed as 9 b \ TtT jd = 
#6 ks © ^l(T CT s)^ implying 



det(^|r CT T d 



det{6' b 



011 



(17) 



In order to compute 0' b \^ Ti7 ^± it should be remarked that the derivative of the equation d b 8 b = X o0 b 
is d b 6' b = X' o9 b 6' b , leading to 6' b (X(a)) = X(8 b (o~)). As the one- dimensional space (T^E)- 1 is spanned 
by X(a), it follows that 



w,, 



bl(T CT S) J 



e:, 



\X(a) 



Consequently 



det(9' b 



X(9 b (a)) 



exp 



\X{9 b {a))\ 
\X(a)\ 



du div(A")(0 u (o-)) 



X( 
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Replacing this in (|T4|) proves ( !T5|) . 

The following notation will be useful 



□ 



dh(cr) 



det(6' b 



f>|T CT £, 



X(B b (a)) 



ex P \ 2 



From ( IT5|) . it follows that the map U defined on T> by 



dudw(X)(9 u (a))) X(a) 



(19) 



(U<f>) b (a) = d b (a) (j)(9 b (a)) , 4> £ V C L (T ) , (20) 

extends to a unitary from L 2 (T d ) to L 2 (IR) ®L 2 (E, z/). The variable 6 is the rescaled energy difference 
w.r.t. the reference quasi-Fermi surface E. Expressing this in terms of W b (see equation (Q), leads 
to (U(p) b (a) = \X(a)\^(W b (f))(a). The inverse, acting on ip e L 2 (R) <g> L 2 (E, u), is given by 



(W»(fc) = d 6 (0_&(fc)) _1 i>b(9-b(k)) 



/(£(*))• 



Note that £Y is unitary. The expression if = UHqU* = U.J- ~HqT*U* will be called the REF represen- 
tation of Hq. Any operator in the REF representation will carry a tilde. The operator (Bip) b = bip b is 
the rescaled energy. Its conjugate operator clearly is A with (Aip) b = -d b ip b . Both of these operators 
are unbounded and have the standard self-adjoint domains. The following result states that these 
notations are consistent with the above. 



Proposition 3 The following relations hold 

UH W = r\B)®lv, Uf{H )U* 



B 



UAU* 



A. 



Proof: The only point to be checked is how the commutation relations of Hq and A, as proved in 
Theorem (TJ are implemented under hi. The first identity results from ( 1201 and 

/- 1 (6)0(^(cr)) = £{d b {a))(j>{d b {a)) = (%A %{<?)) . 



The second formula is obtained from the first one through (unbounded) functional calculus. The 
third one follows from 



(A <g> lU<p) b (a) 



% 



^2 e 



ibA 



AA<j>)(a) = (UA<f>) b (a) 



where IA is expressed in terms of the unitary group W& = e lbA up to the factor |X(cr)|2 which does 
not depend on b. □ 

It is worth comparing the previous construction to the usual one used in scattering theory on 
]R d , where Hq = —A is the Laplacian acting on L 2 (W d ). Then, the (unitary) Fourier transform 
T : L 2 (R. d ) !->■ L 2 (R d ) diagonalizes H , that is, FHqF* is the operator of multiplication by S(k) = k 2 . 
This function has only one critical point at k*_ — corresponding to the minimum of energy E_ = 0. 
The vector field X is defined as in (j3J), now with k e M. d . Let the reference energy be E r = 1 
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so that the (quasi-) Fermi surface S is the unit sphere § d_1 . Furthermore let F(E) = 2E, which 
vanishes at the only critical value. Then X(k) = k and f(E) = f x * = ^ln(E). The flow is 

8b((j) = e b a. As div(X) = d, it follows that db(cr) = e^ db . Therefore the unitary transformation 
U : L 2 (R d ) -»■ L 2 (M) <g> L 2 ^" 1 ) to the REF representation is given by 

(U(f>) b (a) = cf>(e b a) . 
This transformation is discussed and used, e.g., by Jensen [Jen] and also [KR3J. 



2.5 EF representation 

Another natural useful representation is the energy and Fermi surface (EF) representation. A local 
version of this representation is used in the paper by Birman and Yafaev [BY] . It is associated with 
the unitary map V : L 2 (T d ) -»■ L 2 ([£_, E+]) ® L 2 (£, v) defined on V by 



I (let {O'f , - 



\v£{e }{E ^ay)\-. 

The unitarity follows directly from (116)) . It is related to the unitary operator W as follows 



W)b{?) = -^t -J—t mmV) ■ ( 21 ) 

F{E)z \X(a)\2 

The EF representation of an operator on L 2 (T d ) is then obtained by conjugation with V. It will 

carry a circle instead of a tilde, such as H = VH V*, A = VA V* and so on. Any operator that is a 
direct integral in the REF representation is also a direct integral in the EF representation. The first 
example of this type is the Hamiltonian Ho itself: 

(H <f)) E (a) = EM*)- 

More generally, given any fibered operator O = J dbO^'m the REF representation, its EF represen- 

tation is given by O = J dEO E with E = Of/ E \. Another example will be the scattering matrix 
below. The dilation operator in the EF representation can be easily deduced from (T2~T|) : 

(A<p) E (a) = F(E)-d E <l )E (a) + ^F'(E)M*), 

o 

where <fi is in the domain of A, in particular, its derivative is square integrable and vanishes at the 
boundaries of [£L, E + }. 
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2.6 Boundary values of the free resolvent 

Let A C Z d be a finite set. Eventually, A will be the support of the perturbation. Associated with A 
is the subspace £ 2 (A) = C' A L Let II* : C' A ' — > £ 2 (Z d ) be the canonical injection obtained by extending 
elements of £ 2 (Z d ) by zero outside A. It is a partial isometry such that 11*11 is the |A|-dimensional 
projection in £ 2 (Z d ) onto the subspace of elements supported by A, while II II* = 1 C |A|. The finite 
volume Green matrix is defined by: 



This is a matrix of size |A| x |A|. If A = {0} it will be called the Green function. An important basic 
fact about the Green matrix is its Herglotz property, that is, -SmG^z) = i{G^{z) — (z)*)/2 > 
for Qm(z) > 0. This implies, in particular, that G^(z) is invertible for 3m(z) ^ 0. The boundary 
values of G^(z) on the real axis will be analyzed in this section. Gieseker, Knorrer and Trubowitz 
[GKT J studied thoroughly the Fermi surfaces for dimensions d > 2 and for generic periodic potentials. 
They showed that it is an algebraic variety and constructed a compactification. They also investigated 
the nature of the van Hove singularities, which, in two dimension produce a logarithmic divergence 
of the density of states, namely the diagonal elements of 3m G$(E - zO). For d = 2 and the discrete 
Laplacian, these limit behaviors can also be read off the explicit formulas for the Green function 
given in [Ecoj . but for d > 3 only numerical results and toy models seem to be known. 

Proposition 4 Let d > 3 and let £ be analytic. The weak limits G^(E ± zO) = lim^o Gq {E ± ze) 
exist. Furthermore: 

(i) Away from the critical values of S, the map G^(E ±z0) is real analytic. At the critical 
points it is Holder continuous. 

(ii) QmG^(E — zO) = — ^smG^(E + t0) vanishes on (— oo, E-\ U [E+, oo). It is a positive matrix with 
nonzero diagonal entries on (E_,E + ). 

(iii) The map E G M t— >■ ?ReG^(E) is negative and decreasing on (— oo, EJ\ and positive and decreasing 




U(z- Ho)' 1 n* . 





(v) There are matrices N± < such that 
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Proof: The proofs given below are detailed extensions of the work of van Hove [VH]. For m, n G A, 
the matrix elements of G^(z) are given by 



(m\G*(z)\n) = (m\(z - Ho)- 1 ^) = f 

J j° 



d d k e l( - n - m) - k 



l Td (2n) d z - £(k) ' 

(i) Outside the critical values: By construction the matrix G^(z) is holomorphic for z ^ <j(Hq). 
In particular, since the spectrum of Hq is the interval ct(Hq) = [E^,E + ], it follows that the map 
E£l\ [E-, E + ] \-> G$(E) is real analytic and converges to zero at ±00. Moreover, its derivative 
is negative. In particular, if the limit of this matrix exists at E±, this limit is a negative matrix at 
E- and a positive matrix at E + . Now, since S is analytic, it follows that it has only a finite number 
of critical points and it admits a holomorphic continuation in (T + iM.) d in a small neighborhood of 
the form B n = {k + ik G (T + iR) d | maxi<j<d \k%\ < rj}. It follows that, for e > small enough, 
the manifold defined as the set T d = {k + leVS(k) | k G T d } is entirely contained in B v . Using the 
Cauchy formula, it follows that 



(m\G*(z)\n) 



d d k' pi{n—m)-k' 



(2ir) d z-S(k') 



Since k' G T d , it follows that k' = k + ie¥£{k) for some k G T , so that, using a Taylor expansion, 

QmS(k') = e|W(A;)| 2 + C(e 2 ) . 

Consequently, if E G [E-, E + ] \ £(§*) is not a critical value, there is p > such that, if \z — E\ < p, 
the distance of distjz, S(T d )} > does not vanish. In particular, Gq (z) extends as a holomorphic 
function of z from 3m(z) < to a neighborhood of E. In particular, the boundary value G^(E — tO) 
is analytic in E in [E_, E + ] \ S(S*). A similar argument applies to G^(E + tO). 

(ii) Partitioning: For any k* G S*, let Bs(k*) be the open ball centered at k* of radius 5 > 0. 
Let also Bs/2(k*) be the closed ball also centered at k* of radius 6/2. Let U ICS be the open set 
obtained by removing from T d the union of the balls Bs/ 2 (k*), k* G S*. It follows that the family 
{U rcs } U {B$(k*) I k* G S*} is a finite open cover of T d . Let then {x rcg } U {xk* \ k* G S*} be a smooth 
partition of unity associated with this open cover. The previous integral can be decomposed into a 
sum 

r rl d k p t{n—ra)-k 

(m\(z - H )-'\n) = G ICS (z) + G k *(z) , G k ,{z) = ±^ Xk *(k) . (22) 

fc , eS , >>B g (k*) [ Z7T ) Z- t\K) 

The contribution G rog is regular because the integral vanishes around all critical points. Using the 
coarea formula and the results of Appendix it follows that G rcg is holomorphic in the complement 
of the spectrum of H and its boundary values are smooth everywhere on the real line. 

(iii) Non extremal critical points: The boundary values of the Gfc*'s, however, may not be 
smooth because of the contribution of the critical point. Let k* be one of the critical points of 
signature d = {d +1 dJ) with 4^0 and in the following G* = Gk* will denote its contribution to 
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the previous decomposition. If 8 is small enough, the Morse lemma |Nic] implies that there exists a 
neighborhood U of k* containing E>s(k*) and a diffeomorphism ip : Bg(0) — > U such that ip(0) = k* 
and = £ o ip is quadratic: 

i=l j=d + +l 

for \\k\\ < 8 and where E* = £{k*). This diffeomorphism has a Jacobian matrix J = y?'(0) sat- 
isfying Jdiag(ld + , — la ) J* = In particular, the Jacobi determinant of ip stays close to 
| det(£"(/c*))| -1 / 2 over the neighborhood £/ and is a smooth function. It follows that the integral 
defining G> is given by 

/fjdh. p i(n-m)-(f(k) 

It will be convenient to use the following polar variables 



ki = r + u + if 1 < i < d + , fcj = r_w_ if d + < j < d , 

were r± > are the radial variables and u± G § d±_1 the angular ones. It follows that 



GJz) 



ri + 1 dr + r d _ 1 dr_ F(r + ,r_ 



r 2 +r 2_ <S 2 (2n) d z - E* - \{r\ - r 2 _) ' 



where F is a smooth function with support inside the disk r\ + r 2 _ < 8 2 given by 
F(r + ,r_) = / dw + eL;_ |det(^(^)) | Xfc* (<^W) e^~ m ^ 



p(e)de 



Equivalently can be expressed as 
where p is defined by 



z — E* — e 



J r 2 ++r 2_ <5 2 {2n) d V 2 

If e > 0, the usual rule followed by the Dirac distribution 8 leads to 

= / 7Tu r d - l (e + r 2 Y d +- 2 ^ 2 F(V^T^,r) (23) 
Jo l 27r ) 

For e < 0, a similar formula holds by exchanging d + with d- and F(r,r') with F s (r,r') = F(r',r). 
The previous expression shows that, if oL > 2, the Lebesgue dominated convergence theorem implies 
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that the limits p(±0) exist and are equal. In particular, p is continuous at e = 0. Moreover, since 
d > 3, if d+ = 1, then d_ = d-l > 2. Then r^-^e + r 2 )^" 2 )/ 2 = r ^-2( e + r 2)-i/2 < r d-s/2 lowing 
that, again, p is continuous at e = 0. 

Equation (I2"3"j) also shows that p is differentiable for e ^ 0. Moreover, its derivative is given by 
the sum of two terms p\ + p' 2 with 

p[(e) = ^[^y r d -\e + rr +/2 - 2 F{V^Te,r), 

P ' 2 (e) = i£* r d -\e + rr +/2 ~ 3/2 ^F(V^ + i,r). 

The same argument as before shows that, if ci > 3, p' 2 admits a finite limit as ±e 4- 0. However, these 
two limits may not be equal, if F ^ F s . On the other hand, if d > 5, p[ also admits limits and the 
two limits coincide. For d = 3,4, however, it follows that d + < 4 so that p[ may diverge at e — > 0. 
Nevertheless, the integrand can be bounded by 

which is integrable if a > 1/2 for ci = 3 and a > for d = 4. Hence in both cases, there is A' > 
such that 

S | - l*)l < ^e-, 

showing that p is Holder continuous at the critical points. Using the Plemelj-Privalov theorem 
(Lemma [TOl of Appendix |A|) , it follows that the same is true for the boundary values of G*. 

(iv) Near the extrema: The behavior near the maximum or the minimum can be treated similarly 
so that it is enough to consider only the minimum at fc* . Again by the Morse lemma, there is a 
neighborhood U of k*_ containing B$(k*_) and a diffeomorphism (p : Bg(0) — > U with <p(0) = k*_ and 
such that £ o cp(k) = i?_ + (1/2) 5^ i=1 kf. Introducing the polar coordinates r = ||fc|| and u G S d_1 
so that k = ru, the contribution G-{z) = Gk*_{z) is given by the integral 



rS d-l j i(n-m)-kt fP( r ) f 

G-{z) = / — jU, F(r) = / du\ V \ru)\ X ^(ru))e^- m >^- k -K 

with x_ a smooth function with support in U(k*_) which is equal to 1 on the ball ||A; — k*_\\ < S/2. 
In particular, F is smooth and bounded in < r < 5, it vanishes in a neighborhood of r = 5 and 
all its derivatives have a limit at r = 0. Consequently, the integration domain can be extended to 
[0, oo) without change. At this point two remarks should be made: 

(1) F(0) = | de%'(0))| (S^ 1 ! > and the Morse lemma shows that | de%'(0))| = det(£"(£;* )) 1 / 2 . 

(2) The expression e l ( n ~ m )- fc I j s the matrix element |A|(m|M?|n) of the projection matrix M^. 
The change of variable e = r 2 /2 yields 
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Since d > 3, the function e G [0, oo) i— > e^~ l F(y/2e) is continuous and vanishes at e = like e d ^ 2 ~ 1 . 
Hence it can be continued as a Holder continuous function on the entire real line with support in 
[0, y). Consequently, thanks to the Lemma [TUj GL(E±zO) is also continuous w.r.t. E. In particular, 
it has a finite value at E — EL. Since the other contributions to are regular near E_, G^(E ±i0) 
is also a Holder continuous function of E near E = E_. 

All contributions in equation ( |22l) other than G- being analytic near E = E_, it follows that any 
singularity of G^(E ± zO) near E = E_ is coming from G-. In addition, the imaginary part of the 
other contributions to Gq vanishes on the real axis at E = E_ since Hq is selfadjoint. Hence the 
only contribution to its imaginary part is coming from GL. Thanks to the Lemma [TU1 it follows from 
f )24p that this imaginary part is exactly 

SmG?(E ±8 0) = (2(E - E.))^\A\ M n + 0((E-E_)§). 

On the other hand, the real part can be estimated by considering the subdominant contribution of 
G- given by the difference 

GL(E±zO)-GL(E_) = 2(E-E_) Ll^L [°° de SMLllb^L . ( 25 ) 

The same argument as before shows that the integral defines a continuous function of E on the 
real line if ti > 5. Consequently, i? G 1 K G^{E ± zO) is continuously differentiable in a small 
neighborhood of E = E_. Since the derivative is negative outside of the spectrum of Hq, it follows 
that the claim (v) of the Proposition H] holds for d > 5. 

For d = 4, equation (125]) shows that Lemma [TT1 applies. For indeed, the function e G [0,oo| h- >• 
E(y / 2e) G C is smooth because the Taylor expansion of E(r) near the origin contains only even terms 
and E(0) ^ 0. Consequently 

i(n-m)-k*_ jp(r\\ 

G_(E±*0)-GL(E_) = (E-E_)ln(|E-E_|) / ] + 0(|E - E_|) . 

Since all other contributions to the real part of Gq (E ± zO) are regular at E = E_, it follows that 

c\et(f"(k* )Y/ 2 

G?(E ± .0) = G?(E_) + |A|M° 1 ( 2 ~)i x ' & " ln d E " E -\) + ~ " 



At last, for d = 3, returning to the variable r = V2e, equation (|25|) becomes 
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3 i(n-m)-fc* roo F(r) 



i(n—m)-k_ roc 

G_(E ± ,0) - = 2(E - E.) / rfr 

The integral on the r.h.s. can be decomposed into two contributions 



E - E_ - f ± tO 







with 



oo 



h(E) = F(0) / dr - - - . ■ - , I 2 (£) = / dr fflj - . (26) 



o £ - £?_ - £ ± «0 Jo £ - E- - f ± tO 

The first part can be computed explicitly to give 

/ x ttF(O) 

^2{E~E~) 

This part is singular and gives a nontrivial contribution to G^(E ± zO) of the form 

G?(JS ± ,0) = ± . V2(E - E_) |A| ^ |S 1 + / a (£0 , 

where 12(E) comes from the contribution of I2. Since the matrix M 1 } is a projection, the singularity 
does not contribute to the real part of this expression. On the other hand, the integral I 2 can be 
treated by using two remarks: (a) F(r) - F(0) = 0(r 2 ) at r = 0, (b) for < E - E_ < ^ and, 
for r > 5, the integral converges to a smooth function of E. Hence the contribution of the integral 
coming from r > 5 does not produce any singularity, while the contribution for r < 5 is regular at 
r = leading to a contribution that is continuous at E — E_ thanks to the Lemma [TQ1 This finishes 
the proof of Proposition HI □ 



2.7 Localized states in the REF representation 

The REF representation of the localized state at site m G Z d is ip m = UF \m). The states (Vv) m ez d 
form an orthonormal basis in L 2 (R) (g> L 2 (S, v). More explicitly, they are given by 

lM<0 = T~~7d d b {a) e™ 6 ^ , (27) 
(2tt)2 

for any o G S avoiding §. It will be convenient below to consider ip m ^ as a state in L 2 (S, 1/). These 
restricted localized states are not normalized, but their norm is independent of m: 

This norm as well as scalar products between these states are linked to the resolvent. 
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Lemma 2 The following holds 

(ipnA^L^) = (n\ T%m((r 1 (b)±iO-H )- 1 ) \m) . 

Proof: Thanks to the coarea formula and the Plemelj-Privalov theorem (see Lemma [TO]) 

/ Mdo) 1 l(m -n).a 

h E (2tt)« \VS(a)\ 

Always using E = the map 9b '■ £ — > is a diffeomorphism. Thus the associated change of 

variables gives 



{ n\%m((E±iO-H o r 1 )\m) = T vr / |det(^[ TjS )| - W 

=F 7r /" i/(do~ 



V£(9sW)| 

Now the formula follows from the definition of ip m b and ( 1T91) . □ 



Corollary 1 Let us introduce the operator IZb = ^ m eA ( m l mapping £ 2 (A) = C' A ' onto t/ie 

subspace of L 2 (S,z/) spanned by the (Vv>) m eA' ?7ie range o/ zts adjoint TZ\ = Y^meA \ m )(i J m,b\ is 
denoted by Tb C C' A '. Further let IL^ &e a partial isometry from Tb onto the subspace of L 2 (E 
spanned by the {ip m ,b) meA - 
(i) TTie following holds 

Uplb = 3m Gq - zO) , 6 = /(£?) . 

7T 



(ii) i/Pfe denotes the orthogonal projection in C' A ' onto t/ie subspace Tb, th 



en 



(iii) T/ie mop 6 G R i— >• 7?.;, G i3 (C' A ', L 2 (S, z/)) is norm continuous. 

Proof: (i) is a re-phrasing of Lemma |2] and (ii) just the usual polar decomposition, (iii) Since TZb has 
finite rank, the norm continuity follows form the strong continuity. In turns the strong continuity 
follows from the continuity of the inner products (^blVw) l 2 (y,,v)- The latter property follows from 
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Lemma H] and from the continuity of F, f 1 and the imaginary part of the Green function (see 
Proposition 0}, with respect to E or to b. □ 

It is worth remarking that Pf, = n^Ilfe and Ub = HbPb- Furthermore, G^(f~ l (b)±tO) commutes 
with Pfe. The next lemma is a technical result which will be needed to deal with threshold singularities 
in dimension d = 3. 

Lemma 3 Let d > 2 and let the extrema of £ be isotropic in the sense that £"(k±) is a multiple of 
the identity matrix. Then 

lim } </w = e* m - fc i ^ ± , (28) 

where ip± G £ 2 (£, v) are normalized states given by 

if> ± (o) = CexpQ^ du(div(X)(e u (a))±d)\ \X(<r)\? , C>0. 

Proof: Let us first argue that ip± are well-defined and normalizable. The isotropy hypothesis and 
(|HJ) imply div(X)(fc^_ + k) = =f d + 0(\k\). But 9 u (a) converges to k± as u — > ±oo at an exponential 
rate. Therefore the integral in the exponential exists and hence ip± are well-defined. Thanks to the 
definition ff27|) of ^) m .b and to equation (fl9j) the conclusion of the lemma follows. □ 

Remark 1 Without the isotropy assumption that £"(k±) is a multiple of the identity, the l.h.s. of 
(j2"8"j) does not converge to a state in L 2 (E, v). □ 

Lemma [2] also implies that the states ipm,b are in general not orthogonal in L 2 (S,z/), not even 
linearly independent as show the next results. 

Lemma 4 Using the notation of the Corollary HJ one has: 

(i) F b = Ran(3mG?(£ - i0)) whenever E_ < E < E + and b = f(E). 

(ii) In any interval o/R not containing a critical value f(£(S*)), there is a discrete subset without 
accumulation points outside of which the dimension of is constant. 

Proof: The first result follows directly from Lemma [2j Since F does not vanish on (E_,E + ), the 
image is also the image of QmG^(E — zO). In particular, 

dim(JS) = Rank(%mG^(E -tO)) , b = f(E). 

As E i— > ^smG^(E — %0) is real-analytic away from the critical values £(S*), the statement (ii) follows 
from analytic perturbation theory. □ 

The next question concerns whether the rank is indeed changing as a function of E. In addition, it 
is important to have examples leading to a non-maximal typical rank because this will allow produce 
embedded eigenvalues later on. The following result will help to construct such examples. 
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(v\$SmG%(E -iO)\v) = 7T I 

Jt,, 



Lemma 5 With the hypothesis of Lemma HI the orthogonal complement in C' A ' is given by 

Ft = Ker(9fmG£(£7 - zO)) , b = f(E) . 

In addition, a vector v = (f m )meA belongs to J-'jj- (with b = f(E)) if and only if its Fourier transform 
^(&) = ^2meA v m eimk is vanishing identically on the energy surface 

Proof: The first relation comes directly from the first result of LemmaHl As in the proof of Lemma[21 

VE(d(r) \v(o-)\ 2 
(27r) d \V£(a)\ ' 

In particular, since QmG^(E — zO) > 0, it follows that v G Ker(QmG^(E — zO)) if and only if v 
vanishes ^-almost everywhere on E#. But since v is a trigonometric polynomial, it has to vanish 
everywhere on E^. □ 

Proposition 5 If E G [E_,E + ] is not a critical energy, a vector v = (f m ) m( =A belongs to , with 
b = f(E), if and only if there is w G £ 2 (Z d ) such that II* v = (El — H )w. Then w admits a 
Fourier transform having the same degree of regularity as S. More precisely, if £ is of class C r with 
rGNU{oo}U{w}, so is the Fourier transform of w. Moreover, if £ is a trigonometric polynomial, 
then w has a finite support and the previous claim applies to any E G E + ). 

Proof: A direct calculation shows indeed that, if = (El — Hq)w with some w G £ 2 (Z d ), then 
(v \QmG^;(E — t0)\v) = 0. Conversely, if v G J 7 ^, then v is a trigonometric polynomial vanishing on 
He thanks to the Lemma [5l Therefore, since E is not critical, the energy surface is a smooth 
manifold and w(k) = v(k) / (E — £(k)) is analytic in a neighborhood of this surface. Since it is analytic 
outside as well, the result follows. The Fourier coefficients of w defines a vector w G £ 2 (Z d ), which 
decay exponentially fast at infinity by analyticity. In addition, if £ is a polynomial, it follows that w 
is a polynomial as well by Hilbert's Nullstellensatz, which applies even if E is critical. Therefore w 
has finite support. □ 

The last proposition suggests examples of situations for which the kernel is trivial or not. 

Example 1 If A is reduced to one point, then p n (E) being positive for E_ < E < E + by Proposi- 
tion HI the kernel is trivial. □ 

Example 2 Let H be the discrete Laplacian in dimension d = 2 with band function £(ki,k 2 ) = 
008(h) + cos(k 2 ). If A = {(0, 0), (-1, 0), (0, -1), (1, 0), (0, 1), (1, 1), (1, 2), (2, 1)} and w(k) = e^ - 1, 
then neither Hqw nor w are supported in A. Nevertheless (| + Hq)w is supported by A. Hence for 
the energy E = — | = f(b) the set is nonempty. □ 

Example 3 Again Hq is the discrete Laplacian in dimension d = 2. Here we choose the set A = 
{(0,0), (-1,0), (0,-1), (1,1), (1,2), (2,1)}, then w(k) = e lkl e lk2 - 1 is supported in A and so is 
£(k)w(k). Therefore for any E = / _1 (&) the vector v = (E — £)w is the Fourier transform of a 
compactly supported vector. □ 
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If £ is a trigonometric polynomial, more can be said. In such a case there is a finite set S C Z d 
such that £(k) = Ylsei d £s£ lsk with £ s ^ •<=>- s G S. This set S is called the support of H . 
Since H is self-adjoint, £ is real-valued so that £_ s = £ s for all s G S. In particular, 5 is invariant 
under the parity map s h-> — s. The S-interior of a finite set will now be defined as the set of its 
points that cannot jump to the outside using hopping terms from H : 

Definition 1 Let A and S be subsets ofZ d . Then the S-interior A s of A is the set of x G A such 
that x + S C A. 

Example 4 If £ is a trigonometric polynomial with support S and if A is a finite set with nonempty 
S'-interior, then any w supported by A 5 satisfies (E — H )w(x) = outside A. Hence the dimension 
of J 7 ^ is at least equal to the number of points in the S'-interior of A. □ 

Proposition 6 Let A C Z d be finite and let £ be a polynomial. Then, there is a finite subset 
V = V(A) in the spectrum of H , such that dim(J r f ^) is constant for E outside ofV, if f(b) = E. 

The proof of this proposition requires several steps that are described in the next four subsections. 
2.7.1 Prime vectors and convexity 

Lemma 6 (Prime vectors) A vector a = (a±, . . . , a^) G Z d is called prime if it satisfies one of the 
following equivalent definitions: 

(i) Any A > such that Xa G Z d must satisfy A > 1. 

(ii) The greatest common divisor of the coordinates of a is equal to 1 . 

(iii) The map a : x G Z d i— > a ■ x G Z is onto. 

Proof of the equivalence: (i)=>(ii) Let p > 1 be the greatest common divisor of the coordinates 
of a. It follows that a/p G Z d . Therefore 1/p > 1 implying that p = 1. 

(ii) =>- (iii) Let the greatest common divisor of the coordinates of a be equal to one. The map <p a is 
a group homomorphism. In particular, its image is a subgroup of Z. Therefore there is an integer 
p > 1 such that this image coincides with p7L. The coordinates of a are given by a, = a (ej), where 
{ei, . . . , ea} denotes the standard basis of 7L d . Thus, there are integers 6, such that a, = pbi for all 
i's. Since the greatest common divisor of the coordinates is 1, it follows that p = 1. 

(iii) =>(i) Let a be such that 4> a is onto. Then there is y G Z d such that 4> a {y) = 1- Let A > satisfiy 
Xa G Z d . It follows that A G Q and that 4>\ a {x) = X(j) a {x) G Z. In particular, 4>\ a {y) = A G Z, 
showing that A > 1. Hence (i) holds. □ 

Definition 2 A {positive) half-plane in 7L d is a set of the form 'K^m = {x E Z d \ a ■ x > m} where 
a is a prime vector in Z d and m G Z. The associated (oriented) affine hyperplane 0-C am is defined 
similarly as c K a ,m = {x <E Z d \ a ■ x = m} . 
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Definition 3 Let A C Z d . 

(i) A prime vector a is a A-direction if there is m G Z d such that A C J{+ . 

(ii) An oriented affine hyperplane "K a ,m {resp. half-plane m ) is called a contact hyperplane (resp. 
a contact half-plane) for A whenever A C JC^ m and A fl ^K ajm 7^ 0. 

(iii) The convex hull of A, denoted by Conv(A), is the intersection of all its contact half-planea. If 
no prime vector is a A-direction, then Conv(A) = Z d . 

(iv) A is convex whenever it coincides with its convex hull. 

(v) A finite convex set is called a polytope. 

It is easy to check that a polytope has a finite number of contact hyperplanes. 

Lemma 7 Let a be a prime vector in Z d . Then there is a matrix A G SL(d, Z) ; such that <fr a o 
A(n, y) = n for all (n, y) G Z x Z d_1 . 

Proof: Thanks to the Lemma [6], there is a vector bi G Z^ such that <p a {bi) = 1- On the other 
hand, Ker(0 a ) is a subgroup of Z d and therefore it is free. In particular, it admits a basis {62, • • • , bd} 
(minimal set of generators). Clearly the vectors bj can be chosen to be prime. Then let A be the 
dxd matrix with column given by the family {61, . . . , bd}. By construction, A has integer coefficients. 
Moreover, b\ = Ae± and therefore a (Aei) = 1. In addition, <p a (Aej) = 4> a (bj) = whenever j > 2. 
Consequently, if y G Z d is such that Ay = 0, then yi = and J2j=2 vfij — 0. Since Ker(0 a ) is 
free, this implies that yj = for all j's. Hence A is one-to-one. On the other hand, if x G Z d , then 
x — <f) a (x)bi G Ker(0 a ). Therefore there are y 2 , ■ ■ ■ , Vd £ Z such that x — 4> a (x)bi = ^2j- 2 lljbj- Setting 
Ui — 4>a{ x ) and y = (2/1, ■ • • ,yd) it follows that x = Ay. Hence A is also onto. Therefore the matrix 
A is invertible and its inverse has also integer entries. In particular, changing the sign of one of the 
fej's if necessary, det(A) = 1. □ 

2.7.2 The A-Fourier transform 

Let s G Z d and let A G SL(<i, Z). Then let T s and Ua be the operators acting on £ 2 (Z d ) defined by 

(T„V) (*) = V0& - s) , (W) (*) = ^(A-'x) , v e • 

Then both T s and Ua are unitary operators. Moreover, if s, i G Z d , one has T s+t = T s T t and T = 1. In 
a similar way, if A, B G SL(d, Z) then U aUb = ^Tas and t/id = 1. In particular, (t^) -1 = U A -i = U* A . 
In addition, UaTsU^ 1 = T As leading to 

H = ^2^sT s =>- UaHqU a 1 = £-a~ 1 s T s , 

s£S seAS 

Thus L^4 changes the support of H from S into AS 1 . Let now a be a prime vector and let A G SL(d, Z) 
be chosen to satisfy <fi a (Aei) = 1. The partial Fourier transform Ta is the unitary transformation 
from £ 2 (Z d ) into £ 2 (Z) <g> L 2 (T d-1 ) defined by 
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(^)n(p) = Yl H A (n,y))e tp ' v , ^et 2 (Z d ), peT^ 1 . 

It follows that 

W)„(P) = E ^(p)(^AV)n-,(p), ^(p) = E ^)e^. 

re<f> a (S) teZrf- 1 , (r,t)&A- 1 S 

Since S 1 is finite, each of the S^s is a trigonometric polynomial. Moreover, 5* is invariant under 
the reflection s !->■ — s, so that (j) a (S) C [— r a ,rj if r a (S) = max(j) a (S) and both ±r a G 4> a {S). 
Consequently, since (f) a (A(r,t)) = r, it follows that £^ ^ if and only if r £ a (5')- In particular, 

^ o. 

2.7.3 The case of convex A 

From the definition of the S'-interior, it follows immediately that the S-interior of a convex set A 
coincides with the intersection A 5 = f]( a m) m+r a (s) over ^ ne P a i rs ( a , m ) such that 3i+ m is an 
oriented contact half-space for A. This leads to the following result. 

Proposition 7 Let A C 7L d be a polytope and let S be the support of Ho. Then the equation (E — 
Hq)w{x) = for x ^ A admits a solution w with finite support if and only if w is supported by the 
S -interior of A. 

Proof: Let a G 7L d be a prime vector and m G Z be such that 3i a ,m is an oriented contact hyperplane 
of A. Let also A G SL(d, Z) be chosen such that 4> a {Aei) = 1. Then for w G £ 2 (Z d ) let w n (p) denote 
the partial Fourier transform (J r Aw)n(j>)- It follows that, since a (S) C [— r a (S), r a (S)], w satisfies 
(E — H )w(x) = for a (x) < m, namely 

(E Sr,0 ~ Wn-r(p) = , V 71 < 1TI . (29) 

\r\<r a (S) 

In the following r a (S) will be denoted by r a . In particular, £_ Ta ^ 0. Since w is finitely supported, 
there is N G N such that w n = for n < —N. In addition, each component w n (p) is a trigonometric 
polynomial in p. Writing the equation (129]) for n = —N — r a leads to 

£~r a (p)W- N (p) =0 U7_7V = 0. 

Proceeding to write the equation (T2"9l for n = — AT — r a + / for I = 1, . . . ,m — 1 + N + r a , gives, by 
the same argument, w n = for n < m + r a . Hence the support of w is contained in the half-plane 
'H^m+ra- Since this is true for any contact hyperplane < K a ,m, the support of w is contained in the 
S'-interior of A. Conversely, if w is supported by A 5 , it follows that (E — Hq)w is supported by A. □ 
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2.7.4 Conclusion of the proof of Proposition [6] 

Proof of Proposition [6t Since A C Conv(A), it follows that the equation (E — H )w(x) = is 
satisfied for x ^ Conv(A). Thanks to the Proposition [7J it follows that w is supported in the S- 
interior of Conv(A). Let R be the orthogonal projection on the subspace £ 2 (Conv(A) s ) and let P the 
orthogonal projection on £ 2 (A). Then both P and R are finite dimensional. In addition, the previous 
equation is satisfied if and only if w £ Ker((l -P){E- H )R) = Ker(R(E - H )(l -P){E- H )R). 
The matrix A(E) = R(E - H )(l - P)(E - H )R is finite dimensional, acts on £ 2 (Conv(A) s ) and is 
a polynomial in the variable E £ [E_, E + }. Therefore its kernel has a constant dimension away from 
a finite set by analytic perturbation theory □ 



3 Scattering by a finite range perturbation 

This section is dedicated to the scattering of a lattice particle by a finite range perturbation. In the 
first part, the Green matrix of the perturbed Hamiltonian will be investigated, in the second part, 
various formulas will be derived for the wave operators, the scattering matrix and the time delay 
operator. At last, the Levinson theorem will be proved. 

Let A C Z rf be a finite subset, with |A| points. Let IT : £ 2 (Z d ) — > C' A ' be the corresponding partial 
isometry (see the introduction of Section 12. 6j) . The perturbation will be a finite rank selfadjoint 
operator V : £ 2 (Z d ) — > £ 2 (Z d ), supported on A, namely V = n*nVTTTI. Hence V is encoded in the 
|A| x |A| matrix V n = UVU*. For convenience, V n will be assumed to be invertible (this hypothesis 
can be dropped if the kernel is eliminated like in Section |2T7|) . The perturbed Hamiltonian describing 
the scattering is then H = H + V. A typical example for a local perturbation is a potential with 
support A, namely V = J2n€A Vn \ n )( n \ with v n £ M\{0}. 

3.1 Green function 

Let G n (z) = U(z — H)~ l U* be the Green matrix of the perturbed Hamiltonian. As for the unper- 
turbed case, it is also a Herglotz matrix which is invertible for 3m(z) ^ 0. The following formulas 
are well-known. 

Lemma 8 For z £ C \ R, 

G u (z) = (G* 1 ^)- 1 -^)- 1 = (1 - Gjfc^T 1 ^*) > ( 30 ) 
Let the T -matrix be defined by 

T(z) = n*T n (z)n, T n (z) = (l-V u G^(z)y 1 V n . (31) 

Then 
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Proof: The resolvent identity yields 
1 1 



V = 1 — V 



z — Hq z — H z — Hq z — H \ z — Hq J z — H 
Applying II and II* from the left and right respectively gives 

Gg(z) = (l-G%(z)V n )G u (z). 

Now Gq (z) is Herglotz and thus invertible since z ^ R. Hence, G 1 (z)^ 1 — V u is also Herglotz and 
invertible, leading to the invertibility of 1 - G^{z)V n = G 1 (z)(G J (z)- 1 - V n ). To prove (J32D, the 
resolvent identity gives the factor 1 — Z _ H V. This operator is invertible, because it is a finite rank 
perturbation of 1 and any element in its kernel is an eigenvector of H — Hq+V with eigenvalue z ^ R, 
namely the kernel is trivial. Using the identity (1 — A)~ X A = (1 — A)^ X {A — 1 + 1) = (1 — A)^ 1 — 1, 
the inverse can be written as 

v\ = i + — - — n* v u n (i — v 



Since 



it follows that 



z — Hq J z — Hq \ z — H( 

n (i-j^jf/) 1 = (i-G^v^y'u, 

When combined with the resolvent identity this completes the proof. □ 
3.2 Spectral analysis 

Because the perturbation has finite range, the essential spectrum of H is given by the essential 
spectrum of Hq. However, H may have some discrete spectrum, which, since H is selfadjoint, is 
given by the simple poles of the resolvent on the real axis. Thanks to Proposition HI it follows 
from equation fl30|) that the only way to get a polar singularity in the Green matrix of H is for 
1 — G J (z) V n = ((V n ) -1 — G J (z)) V n to have a nontrivial kernel (recall that we restrict ourselves 
to the case of invertible V n ). This can be analyzed using the determinant of 1 — G^{z) V u which is 
also called the perturbation determinant |Yafj . Furthermore, if E is an eigenvalue of H, 

multiplicity of E = dimKer ((^ n ) _1 — G J (E ± zO)) . (33) 
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If E (jL [E_,E + ], it is called an isolated eigenvalue while, if E G (E_,E + ), it is called an embedded 
eigenvalue. For E = E±, a non-trivial kernel of (V 11 ) -1 — G J (E±) leads to a threshold singularity 
which will be dealt with below. With any E G R is associated the subspace of C' A ' 

S E = Ker( (K n ) _1 - G J (E) ) . 

Then the multiplicity of the eigenvalue E G M\ [.EL, EL] of i7 is also equal to dirties). The embedded 
eigenvalues are characterized in the next result (where the space Te is the space for E = f(b) 
used in Corollary [Q . 

Proposition 8 Let V have finite range. Then E G (EL, EL) is an embedded eigenvalue if and only 
if Te H Se is non-trivial and the dimension of this intersection is equal to the multiplicity of E. 
If £ is analytic, the associated eigenvectors are decaying exponentially fast at infinity. If S is a 
trigonometric polynomial, the the eigenvectors have a finite support. 

Remark 2 The last statement, namely that the eigenvectors have compact support, was proved in 
a slightly different context in |KVj . □ 

Proof: Let E = G {E-, E + ) be an embedded eigenvalue and v G C' A ' be the associated vector 

in the kernel of (V 11 )- 1 - G$(E ± i0). Because 5sm(G$(E - i0)) > and V n is self-adjoint, this is 
equivalent to having 3m Gq(E — iO)v = and (V n ) _1 f = 9fte G^(E)v, or alternatively v G J 7 ^ C\Se- 
As shown in Lemma and in Proposition for any vector v G J 7 ^ fl Se there is w G £ 2 (Z d ) such 
that v = (El — H )w. Moreover w decays exponentially fast at infinity and has finite support if £ 
is a trigonometric polynomial. □ 

Example 5 Let A be such that J=^ ^ {0} (see Examples M to gj) . Chosing V u = (KeG^(E))- 1 
whenever the latter is invertible leads to a Hamiltonian H with an embedded eigenvalue of multiplicity 
dim(J^). □ 

Example 6 Let us present another way to construct Hamiltonians with embedded eigenvalues, again 
by perturbing a periodic H with a polynomial energy band £ . Let P = 11*11 and Q — 1 — ITl! and 
then set V = —PHqQ — QHqP. Clearly V has finite range. Now H = Hq + V splits into a direct 
sum of PHqP and QHqQ. The former has finite rank and admits a spectrum of eigenvalues with 
eigenvectors supported in A. By the minimax principle, it follows that all these eigenvalues belong to 
[E_, E + ]. On the other hand, QHqQ is a finite rank perturbation of H , so it has the same essential 
spectrum. Hence the eigenvalues of PH P are embedded indeed. □ 

It is natural to address the question of whether embedded eigenvalues exist if V is a multiplication 
operator with support in A. The following result gives a negative answer. 

Proposition 9 Let Hq be a trigonometric polynomial and V = Y2n£A v n \ n )( n \ * s a potential with 
entries v n ^ for n G A. Then H + V has no embedded eigenvalues. 
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Proof: Clearly V u is invertible in C' A l Thanks to Propositions [5] and [BJ if E is an embedded eigen- 
value, then there is a vector v G C' A ' such that (V n ) _1, u — G^(E))v = and n*t> = (E — H )w for 
some w G £ 2 (Z d ). It then follows from Proposition [7] that w is supported in the S-interior A of A, if 
S is the support of H . If A 5 is empty, then v = and there is no embedded eigenvalue. Now let A s 
be non-empty. Because w is supported by A 5 , one has Gq (E)U(E 1 — H )w = Uw so that replacing 
v = U(E 1 — Hq)w shows U(E 1 — Hq)w = UVw and thus (El — Hq)w = Vw. Now the r.h.s is 
supported in A s since V is a multiplication operator. Applying again Proposition [7J shows that w 
is supported in the S'-interior (A 5 ) 5 of A 5 . Now this procedure can be iterated. As the iterated 
S'-interior of A is empty after a finite number of steps, one concludes that w = 0. □ 

Let us next investigate threshold singularities in more detail, now considering only the case of 
dimension d > 3. They appear at either one of the band edges E± whenever (V n ) -1 — G^(E±±tO) has 
a non-trivial kernel, which is equivalent to Se± being non-trivial. Hence dim(iSg; ± ) is the multiplicity 
of the threshold singularity at E±. Again it is easy to produce such singularities by an adequate 
choice of V. Such a singularity can either be a threshold eigenvalue or a threshold resonance (the 
latter is also called a half-bound state) depending upon whether the equation Hip = E±ip has a 
square integrable solution ip or not |Newll UK] . In order to analyze the threshold singularities let 
the following space be defined 

T± = { v G C |A| | v has a zero of order at least 5 — d at k* ± } . (34) 

The definition of 7± is somewhat similar to J 7 ^ as given in Proposition [51 However, 7± is likely to 
be larger than the limit of J 7 ^ as b — > ±oo because it tests zeros only at a single point. In addition, 
7~± coincides with C' A ' for d > 5. 

Proposition 10 Let d > 3 and V be of finite range. Then the multiplicity of E± as threshold 
eigenvalue of H is equal to dimf^Sg^ Pi 7±) . Moreover, the multiplicity of E± as threshold resonance 
of H is equal to dim^^) — dim^^ fl 7±) . In particular, for d > 5 all threshold singularities lead 
to threshold eigenvalues. 

Proof: Let v G S E± . Then V n G$(E±)v = v and V(E ± -H Q )- x U*v = Il*v. Hence there is a threshold 
eigenvector whenever the equation U*v = (E± — H )ip admits a solution ip G £ 2 (Z d ). If so, then 
Hip = E±ip indeed. Using the Fourier transform, the equation leads to the following solution ip(k) = 
v(k)/(E±-E(k)). The denominator satisfies E±-B(k) = ((k - k±) \£"{k±) \ [k - k* ± )) + 0((k- k±) 3 ). 
This quadratic singularity is integrable in dimension d > 3 so that indeed ip G £°°(Z d ) is well-defined. 
If d > 5, then the singularity is also square integrable so that ip G £ 2 (Z d ) is an eigenvector. In 
dimension d = 4, ip is square integrable only if v has a zero at k±. In dimension d = 3 the zero has 
to be of order 2. Combining this leads to the definition fl34l) and to the conclusion above. □ 

3.3 The wave operator as an integral operator 

The potential being finite rank, the Kato-Rosenblum theorem for trace class scattering theory [RS 
lYafj implies that the wave operators 
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Ob = s-lim e lHt e' lHot . 

t— !>±00 



exist and are complete, that is, Ran(fi + ) = Ran(fi_) = P ac (if) where P ac (if) is the projection on 
the absolutely continuous subspace of H . Then the wave operators are partial isometries satisfying 

Q* ± Q ± = 1, Q ± Q* ± = P ac (P) = 1 - P PP {H), (35) 

where P PP (H) is the projection on the pure-point spectrum of H and the last equality holds because 
there is no singular continuous spectrum. In addition, if one sets Q(t) = e lHt e - lH ot^ then Q(t)e lH ° s = 
e lHs Q(t — s) for all sGi Passing to the limit t — > ±oo yields the intertwining relation 



V±g(H ) = g(H)n ±} g e C 0[ 

Birman's invariance principle [RSt lYafj can now be expressed as follows. The function / : (E_, E + ) — > 
R defined in ([7j) is smooth and has positive derivative. It is therefore admissible for the invariance 
principle so that 

tt± = s-lim P ac (#) e lf[H)t e - lf{Ho)t . (36) 

t— s>±oo 

However, H may have some spectrum outside (EL, E + ) and possibly eigenvalues at E± so that f(H) 
may not be well defined. However, using the completeness of the wave operators, P ac (if) can be 
inserted discounting all eigenvalues. Since H has the same essential spectrum as H , this eliminates 
all ambiguity in the definition. It will allow to derive an explicit formula for Q± = J z Vt±J r * which 
will serve as a tool to calculate the wave operator and the scattering operator. 

Proposition 11 The following formula holds 



p JUt p t(k-n—k'-m) 



JL.I ^ilk-n—k'-m) 
n,m£A 

Proof: It follows from DuHamel's formula and a Tauberian lemma IRSI that 



n± = 1 ± z s-lim / ds e ±tHs V e TlH ° s = 1 ± i s-lim / ds e~ es e ±lHs V e TlH ° s . 

Hence 

((h±-l)(f))(k) = iilim f dse~ es (F e ±lHs V e TlH ° s 7* <p) (k) . 
In the following, the notation VJ m = (l\V\m) will be used. In addition, 
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{m\e^ Hos r\<P) = J -^y e ~ lk '' m e T^£(fc ' )s 0(A ; , ) 
Consequently the previous formula leads to 

((Ot-l^W) = ±i V V hm lim r^e"" (j-e^lO) (*) / ^ e^'™ e^ fc ') s . 
V / /mgA 40 J J Jd [271) 

The integral over s can be performed to give 

((SU -!)*)(*) = lim £ X^,.-"- (jr ^J— ^ H>) WW . 

In the previous expression, it becomes possible to compute the part in the parenthesis. For indeed, 
using the resolvent identity as in Lemma [8] yields 

1 ir= i ir 1 



zl-H zl-Ho 1 - V IL G 1 (z) ' 

and remarking that I G A. Hence, passing to the Fourier space leads to 

('rhf i'>) <*> = S ^^"(-Ki-^gjw)- 1 w. 

v ' neA v y 

Replacing this in the above expression for Q± — 1 completes the proof. □ 



3.4 The wave operator in the REF representation 

In this section, the REF representation will be used to calculate the wave operator Q± = UVt±lA* in 
dimension d > 3. It is an operator on L 2 (R) £g> L 2 (S,z/). From Proposition [TTl the definition (|20|) . 
the change of variables formula f|T5|) and the definition ff27j) of the states it follows that 

«n ± - = lim / m £ I'M ( y/r' ( !?-f-'m ' 

J n,meA J ^ ' ~ l ~ ^ ^ ' 

where stands for the inner product in the Hilbert space L 2 (S,i/) and the integral of b' 

carries over R. Thanks to Corollary [TJ the sums over n and m can be computed to give 

(fe - i Wl = ita/ f n ; I^GfCT' + ^ • (37) 

where 0± = fdb 0±f, with 
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0±, b = lim b 1 b T ^ (f- 1 (b)T^e) |9m G^f' 1 (b))\* U b . (38) 

It is part of the proof of the following result to show that the limit in f|38|) exists and that the 
expression f l3T]) for the T-matrix can be replaced to give 



± , b = b 1 b ((V^-UeG^r^b)) T^\QmG^(f-\b))\y 1 ^mG^r'ib^U,,. (39) 

Theorem 2 Let d > 3 and let V have finite support. In addition, F is chosen as in ([3]) and the 
following will be assumed: 

(i) If d = 3, the threshold singularities have multiplicity at most 1 and any vector w± in the kernel 
°f (y n ) _1 — Cq (E±) has a Fourier transform satisfying w±(k±) ^ 0. 

(ii) If d = 4, there are no threshold singularities. 

(hi) The embedded eigenvalues lie neither on the critical values of £ nor in the set V described in 
Proposition O The corresponding zeros ofb G R i— > (V n ) _1 — Gq (/ _1 (6)) are of first order in the 
real part. 

Then the operators 0±^ are well-defined, continuous in b and uniformly bounded. The wave operators 
are given by 

n± = i+J2^*B^ mG o(r\B))\h^ (±i + ewA - Kt \ d ± . (40) 

Formula ( I40p shows that the wave operator can be calculated in terms of 0± and the dilation 

operator A. It is similar to those obtained by Kellendonk and Richard for continuous scattering 

systems [KR2l IKR31 IKR4] , however, we stress that we also allow for embedded eigenvalues, a fact that 

is closely linked to proving that the inverse in (1391) exists. In fact, if the kernels of (V 11 )^ 1 — 9fte G^(E) 

and of ^smG^(E) have a non-trivial intersection, the associated pole in ( 139]) is attained on the 

i 

orthogonal complement of the range of SmGj(/ _1 (6)) 2 and thus one can prove that the inverse 
on J-;, exists (in the sense of Lemma fT2l) . Energies at which this happens are exactly the embedded 
eigenvalues by Proposition El The intersection of the kernels is supposed to be a regular singular 
point and this allows to argue for the continuity of 0±$. Hypothesis (iii) holds generically for V n 
(within the non-generic situations of embedded eigenvalues). Another comment is that the condition 
(i) imposed for d = 3 implies that the threshold singularity is a threshold resonance. This is because, 
by Proposition [TTJ], w±(k±) ^ implies Se± fl 7± = {0}. Again this is the generic behavior in 
dimension d = 3. It is possible to treat threshold singularities for d — 4, but this is technically more 
involved and not carried out here. 

Proof of Theorem [2J Let us first suppose that 0± are well-defined and bounded with fibers 0±$ 
depending continuously on 6, and then show how (T40|) follows from (!37j) . Thanks to the formulas ([7]), 
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S(9 b (a)) = / _1 (6) = E r + Atanh(fr) and F(f- l (b)) = Acosh~ 2 (6), a bit of algebra now leads to 



m ± - mt = n; la^tr-p))!* / g sinh(t ,l 6)Tt0 («* + «-*) io ± , 



V 



In the previous formula, 0±<p is a vector in the Hilbert space L 2 (M)®C' A L As previously let A = —idb 
be the generator of the translation group in L 2 (M) £g> L 2 (S, v) as well as L 2 (M) ® C' A L Changing the 
integration variable 6' to u = b' — b leads to (0±4>) u+ b = (e tAu O±(j))b- Hence 

((fU-l)0) = J2 n b \QmG*(f- 1 (b))\te*i [ - . -±— ^(e^d ± A . 
V /& i 1 J ii smh(u)=FzO V /& 



Now (140p is obtained from the following identity: 

du 1 u r e 71 ^ /t? 

e 2 e = ±1 + — 

27r sinh (-u) =f iO e nA + k i 

It remains to show the above mentioned properties of the operators 0±j, defined in (1391) . We first 
check that for every b G R they are well-defined and continuous in b. This analytical issue is tied to 
embedded eigenvalues. In fact, away from them there would be nothing to prove because the inverse 
in fl39|) exists by fl33|) and the fact that A + iB is invertible for any invertible operator A = A* and 
non- negative operator B > (see the proof of Lemma [T2l) . Now focussing on embedded eigenvalues, 
let us set 

A b = (V^-neG^f-^b)) , B b = \%mG%(f-\b))\ . 

These matrices A b and B b have nothing to do with the dilation operator and the rescaled energy 
operator, and only appear again in the following lines and Appendix [B] Then A b is self-adjoint and 
Bb is non- negative, and by Proposition H](i) both are real analytic in b as long as is not a 

critical value of S. Now the properties of hypothesis (iii) of Theorem [2] guarantee that Lemma [TBI of 
Appendix [B] can be applied because, in particular, the zeros of b e M i— > A b + iB b are of first order in 
Ab. This lemma implies that 0±b is even analytic for b = f^ 1 {E) away from the set £(§*) of critical 
energies and away from the exceptional set V of Proposition [6j Continuity at the latter points follows 
again from Proposition H] because there are no embedded eigenvalues there. In particular, let us also 
note that the dimension of Ub changes at points b corresponding to energies in V, but in (139]) this 
does not lead to discontinuities due to the factor l^sm Co directly following U b . 

Now we check that the operator 0± seen as a linear map from L 2 (R) £g> L 2 (T,, v) to L 2 (IR) <g> C' A ' 
is actually bounded. This means that we have to analyze the limits b — > ±oo of 0±, which depend 
on the behavior at the thresholds. Here the factor ea + e~2 introduced in fl38l) will turn out to be 
crucial. We start by expanding the inverse in (1391 around the band edge E± using items (iv) and 
(v) of Proposition HI 

(V 11 )- 1 -G^if-W^iO) = (V 11 )- 1 - G$(E±) ti D ± Mle-\ b \ {d - 2) + 0(e-\ b \\\b\e- m ) . (41) 
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If there is no threshold singularity, then by definition (V u ) 1 — G^(E±) is invertible and thus the 
inverse in fl39|) remains bounded as b — > ±00 and the other factors lead to lim^-too 0±j, = 0. If there 
is a threshold singularity of multiplicity dim(iS>E ± ) = 1 in d — 3, the operator (V n ) _1 — Gq (_E±) has 
a kernel of dimension 1 spanned by some vector w±. If, in addition, w± is not orthogonal to the 
one-dimensional projection M±, namely if (v±\w±) = \A\~2w(k±) 7^ 0, then ((V 11 )' 1 — Gq (/ _1 (6) =F 

tO))- 1 = O(e^). As \%mG$(f- l (b))\^ = 0(e~^) (because d= 3) the added prefactor (el + 

assures that O^^ remains bounded as b — > ±00. (Its convergence is analyzed in Proposition [TBI 
below.) 

For d > 5, the term of order e~' b '^~ 2 ^ in (jJTJ is dominated by the terms of order e~ 2 ' 6 '. Hence, 
(|4ip becomes rather 

(V 11 )^ 1 — G^(f~ 1 (b) =f i0) = (V 11 )^ 1 — G^(E±) — N± e~ 2 ' 6 ' + C(e~ 3 ' b ') , 

which follows also from Proposition |H In this case, N± is invertible so that for an arbitrary thresh- 
old singularity ((V 11 )^ 1 — Gg(/ _1 (6) =F z0)) 1 = 0(e 2 ^), even in the most singular case where 
(V 11 )- 1 = G${E±). This is compensated in (J3HD by the other factor (e? +e~?)~ 1 \$sm G" (/ _1 (6))|^ = 
0(e~^ ( - d ~ 1 ^). Consequently 0± is bounded for d > 5. This concludes the proof. □ 

3.5 The scattering operator 

Whenever the wave operators are complete, the scattering operator is defined by: 

s = n* + n_. 

It is unitary and satisfies [S, H ] = 0. Hence, in the REF representation, [S, B] = and thus 
S — J db Sb with unitary operators Sf, on L 2 (£,z/). The intertwining relation and the invariance 
principle fl36l) imply that for any admissible function / with /' > 0, one has 

s-lim e ltf{Ho) n ± e - ltf{Ho) = Q1Q± = 1, s-lim e ltf{Ho) Q± e - rtf{Ho) = Q*tt±. 

The second expression is either S or S*. Let now / be chosen as in ([7]). In the REF representation, 
Proposition [3J then leads to 

s-lim e ttS tt±e- ttS = 1, s-lim e ltS fi_ e~ ltS = S, s-lim e ttS tt + e~ ltS = S* . (42) 

t— >±oo t— >oo t— >— 00 

Using the explicit formula for Q_ given in Theorem [2] now leads to an explicit expression for the 
on-shell scattering matrix. The structure of such formulas (in particular, the EF representation of the 
formula (jHJ) in the proof below) is well-known and has appeared in various guises (see |Newll lYaf] 
for a list of references). 
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Theorem 3 Let the assumptions of Theorem [2] hold. Then the on-shell scattering matrix S b is a 
unitary operator on L 2 (S, v) depending continuously on b and given by 

s b = (i-n*n 6 ) + n*(a-*) (a + 

where the selfajoint L x L matrix C b : PfeC' A ' — > P&C' A ' is defined by 

C b = P^mG^f-^b))^^ 11 )- 1 -MeG*(f-\b)))\ZmG*(f-\b))\^ P b . 

Proof: For any function g the following formula holds e ttB g(A)e~ ltB = g(A — t). The limits t — >■ ±00 
can be taken whenever g has limits at infinity. The function appearing in ( 140]) is of that type. The 
middle formula in equation ( l4~2j) and the expression of Q- given in Theorem [2] leads to the calculation 
of S, namely 

S b = 1+ ^n* ^mG^r 1 ^))^ e 4 (-2) d ±tb , (43) 

K=±l 

Because Theorem [2] states that 0± )b is continuous in b, this formula already shows that S b is contin- 
uous in b. Using equation f l39l) . it now follows that 

S b = l-2tW b ISmG^CT 1 ^))^ ((V u )- 1 -G^(r 1 (b)+tO)y 1 \^mG^(f- 1 (b))\^U b . 
After simplification, one gets 

S b = l-2iW b (C b + i)- l Ii b . (44) 
This allows to prove the claim. □ 

Similar formulas hold for the EF-representation of the scattering matrix. The comments made 
in Section [231 and the results of Theorem [3] lead to (with Ce = C b for b = f(E)), 

Se = S m = (I-W^+W^Ce-i^Ce + i)- 1 ^. 
It is now possible to get results on the asymptotics of the scattering matrix. 

Proposition 12 Let the assumptions of Theorem [2] hold. 

(i) If there are no threshold singularities or if d> 5, then lim^-too S b = 1. 

(ii) If for d = 3 there is a threshold singularity of multiplicity 1 at E± and the extremum of £ at E± 
is isotropic in the sense that £"(k±) is a multiple of the identity, then 

lim S b = 1-2|V±)(V±|, 

b— >±oo 

where ip± G L 2 (E, v) are the states given in Lemma [3l 
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Proof, (i) If there are no threshold singularities, then hin^^oo 0±j, = as was shown in Section |3~H 
As l^rnGg 2 e K 2 is bounded, it follows from (1131) that lim^ioo Sb — 1- F° r d > 5, it has 

been shown that 0±.& is uniformly bounded even in the presence of threshold singularities. As the 

- b 

factor 3m Gq | 2 e K 2 vanishes in the limits 6 — > ±oo, the same conclusion holds thanks to 

equation (1151) . 

(ii) For d = 3, it is assumed that there is a threshold singularity of multiplicity 1 and the extrema 
are isotropic. Starting from ( 1H|) the inverse of + i can be expanded using Proposition H] to give 

(C b + i)- 1 = D ± e^ bl M2({V u )- 1 -^eG u (E ± )+iD ± e- w M^y 1 + 0(e~^) 

where, in the second equality, Lemma M stated below is used. It can indeed be applied thanks to the 
hypothesis stated in Theorem El On the other hand 1ZbV± = ^o,b + 0(e Tb ) so that it follows that 
n^M^n b = ll^cftll -2 \i J o,b){4'o,b\ + 0(e Tb ). Hence Lemma [3J allows to conclude the proof. □ 

Lemma 9 Let P 6 Mat (nxn, C) be a one- dimensional orthogonal projection and A = A* G Mat(nx 
n, C) have a one- dimensional kernel not lying in the kernel of P. Then 

XPiA + iXP^P = -iP + 0(\). 

Proof: This follows from a short calculation using Cramer's rule. □ 

3.6 The contributions of the threshold singularities 

Just as the scattering operator is obtained as a rescaled energy boost of the wave operator in (112]) . 
it is natural to study the dilation operator boost of the wave operator. As the scattering operator is 
obtained as boost of it is sufficient to consider 

R± = s-lim e ltI n.e- ltA . 

t-^±oo 

It ought to be remarked that the limits in R± approach the critical values E± respectively. This is 
because the identity e ttA g(B)e~ dA = g[B + t) holds for any function g. From the definition it follows 
that [R±, A] = so that R± = J daR±^ a with operator fibers i?±,o acting on L 2 (£,z/). Since the 
operator A has continuous spectrum, it follows that lim^-too e ltA Ke~ ltA = for any compact operator 
K. In particular, since f2_ is unitary modulo a compact operator, R± is unitary. Consequently each 
-R±, a is unitary. The following operators are associated to R± in a similar manner as the time delay 
is associated to the scattering operator:, 

T ± = ±-(R ± )- 1 [B,R ± }. 
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Proposition 13 Let the assumptions of Theorem [2] hold. 

(i) If there are no threshold singularities or if d > 5, then R± = 1 and T± = 0. 

(ii) Ifd = 3, if there is a threshold resonance at E± of multiplicity 1 and if the extrema of £ are 
isotropic, then 

R±,a = (l - |^ ± W±|) + JrjH Tr(T ± ) = 1. 

Proof: Statement (i) is a consequence of the proof of Theorem [2] and the arguments below, so let us 
focus on the case d = 3. Equation ( 14"CT|) gives an explicit expression of the wave operator. Each term 
in the sum over k = ±1 in its r.h.s. is the product of three terms that can be treated separately 
under the boost action. The first factor is 

lim n; \QmG*(f-\b))\* e«* = (Z>±)* S Kt± |V±}^±i, (45) 

where S Kt ±x is the Kronecker delta. Details for the proof of the second equality are similar to the 
proof of Proposition [T2"l Taking into acount the Ket {v±\ in ( |45p permits to treat the last factor 
using Lemma [91 This gives 

i 

lim e lt ~ A {vl\ 5 ± e~ lt ~ A = lim (v*\ ( (V 11 )- 1 - G*(f-\b) + lO) Y* l SmG ? (/"'p))! 2 n& ^ 

t— >±oo ~ 6->-±oo V / g5 4- e — 2 

The limit on the r.h.s. is given by i(D±)~^ (ip±\. Replacing these factors leads to the formula for 
R± A - Using the integral over a, the trace of T± is given, up to the sign, by the rotation number of 
the map a6R4 (e™ =f t)/{e na ± %) e S 1 . The latter is equal to ± \. The sign in the definition of 
T± compensate the previous sign leading to the final result. □ 



3.7 The time delay operator 

The time delay operator T is the derivative of the scattering matrix w.r.t. the energy (the notation 
T should not be confused the T-matrix). More formally, it is defined by T = - S] whenever 

S is differentiable w.r.t. to the dilation A. In the REF it becomes 

f = fdbT b: f b = ± 



while in the EF representations it is given by 



E 4 



1 



T= dET El T E = -(S E )- l d E S E . 



i 



The total time delay is the trace of T. The formula given in the following result is sometimes called 
the spectral property of the time-delay [TO| INewl] 
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Theorem 4 Let the assumptions of Theorem [21 In addition, suppose that Te = C A ' for almost all 
E and that there are no threshold eigenvalues. Then, for almost all E e [E_,E + ], 

TrL2{s,v)(T E ) = lim 2 Tr^ (z<i) (- - — j . (46) 

40 \E — xe — H E — xe — H ) 

Remark 3 The condition T E = C A ' for almost all E, implies that H has no embedded eigenvalues 
(see Proposition E]). If H has embedded eigenvalues at energy E with multiplicity n(E), then the 
r.h.s. of equation ( |4"6]) must be modified by subtracting n(E)/e to the trace to compensate for the 
singularity occurring at this energy. One should also be able to deal with a threshold singularity by 
subtracting the adequate contribution. However, no details are provided here. □ 

Remark 4 The previous formula for the total time delay is well-known for potential scattering 
when H is the Laplacian in R d (see |CN] for a list of references). It can be proved by a direct 
calculation in the REF representation (following the lines of |TOj ) or by a computation inspired by 
the Birman-Krein formula |Yaf] for the scattering phase, an approach used below. □ 

Proof of Theorem [4} The following notation will be used IT^ = U.fi E ), Ce = C/(s) etc- From 
equation OH]), it follows that 



l-d E S E = -dEll^CE + tr'llE + Il^CE + tr'dECEiCE + ir'llE-IlEiCE + tr'dEllE 

1% 

The equation IT^II^ = 1 C |A| implies Oe^e^e = ~^e9e^*e- Hence, 

Tr L 2 { x >v) {f E ) = ±Trr B (n E (S E )*U E U E d E S E U E \ = 2 Tr c , A| ((C| + l)- 1 ^^) . 



This can be rewritten as 



Tr L 2 (s } (T E ) = - d E In det ( — - - ) 
i \C E + ij 



Therefore dividing out the imaginary part of the Green function appearing in the definition of C E 
(see Theorem [3]) gives 

Tr L2( E,,)(T E ) = \d E \n det(p s ((\/ n )- 1 - G$(E - zO))P i; (P £ ((V rn )~ 1 - G$(E + zO))P B )" 1 ) 
= 2Qmd E In det(P E {{V u )~ 1 - G%(E - iO))P E ) ■ 
On the other hand, using ( 1321) and the cyclicity of the trace, leads to 

Tv Pm ({z-H)^-{z-H )-') = ^(iLiz-Hor^dV^-G^z))- 1 ). 
Since d z G$(z) = -U(z - H )- 2 U*, it follows that 

Tr e2{zd) ((z - H)- 1 - (z - Ho)' 1 ) = d z Tr c , A| (h^V 11 )- 1 - . 
This leads to the identity. □ 
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3.8 A Levinson-type theorem 



Theorem 5 Let the assumptions of Theorem [2] hold. Further let N = Tr(P pp ) be the number of 
bound states of H , including embedded eigenvalues and threshold eigenvalues. Then 

1 rp ( rri\ AT f - § dim (S^J - f dim (<S E _ ) if d = 3 , . , 

-Tr w (T) + N = | + Q . frf ^ 4; (47) 

where dim(SE±) G {0, 1} is the multiplicity of the threshold resonance. 

Two proofs will be provided. The first one will require stronger hypothesis to show how the 
argument principle combined with the spectral property of the time delay can be used as in most 
standard references \RS\ INewlj . It might be possible to lift these hypothesis with more technical 
effort (see Section [3T91 where this is done for a point interaction). The second proof requires only the 
stated hypothesis. It is based on the approach proposed by Kellendonk and Richard [KRlj and uses 
on the index theorem for C*-algebras. 

First proof of Theorem [5j It will be assumed that there are no embedded eigenvalues and no 
threshold singularities and that T E = C A ' for almost all E (these assumptions hold for the case of a 
single site impurity). The number N of eigenvalues is obtained by counting the poles of the resolvent 
using the Cauchy formula and a contour integration. The contour is given by two circles, one large 
counterclockwise oriented circle 7 around the spectrum of H and a second small clockwise oriented 
circle Y around the spectrum of Hq (but not touching it). Then 

N = I Tr, 2(zd) ((* - H)- 1 -(z- Ho)- 1 ) . (48) 

The resolvent identity implies that the contribution of 7 vanishes in the limit where its radius goes 
to infinity. Then let T converge to the concatenation of the two intervals [E_ + iO, E + + 2O] and 
[E- — tO,E + — z0]. Since it has been assumed that there is no threshold singularity, the regularity 
of the Green function at the band edges implies that the small circle connecting these contours near 
the band edges have a vanishing contribution in the contour integral. Thus 



dE m / 1 1 1 
N = I — Tr £ 2 {zd) I - : - - : — - — + 



l E _ 2m ( >\E + iO-H E + zO-H Q E - 1O - H E-i0-H 0/ 
The formula for the total time delay, proved in Theorem HI gives 



N = / + dETr L ^ v) {T E ) 

■ZK J E _ 



The r.h.s. is nothing but the trace of T expressed in the EF representation, leading to the result. □ 

Second proof of Theorem i As a preamble let us construct an extension the algebra Co(M) of 
the continuous functions on M vanishing at +00 and —00. Then set Coo(M) = C(M) where R = 
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{ — oo}UlRU{+oo} which are the continuous functions having limits at +00 and —00. The evaluation 
map ev is the *-homomorphism defined by ev : g G C(M) h-> ev(g) = (g(+oo),g(—oo)) G C 2 = C©C. 
The kernel of this map is precisely C (K) leading to a short exact sequence —> C (R) Coo(R) — > 
C © C — 7- 0. Next let us consider a two-dimensional version of this extension. Let Coo(IR 2 ) be 
those continuous functions on IR 2 having a continuous limit function built from the limits of / in 
the directions of M 2 . These directions can described by S 1 so that one obtains an exact sequence is 
— > C {R 2 ) Coo(lR 2 ) C(S X ) — > 0. On the other hand the limit points at infinity can also be 
seen as a square the corners of which being given by points of coordinates (±00, ±00). Equivalently, 
Coo(lR 2 ) can be seen as the subalgebra of the C*-tensor product Coo(IR) © Coo(R) generated by 
functions that are continuous at the four corners of the square at infinity (let us point out that 
there is a unique C*-norm on the tensor product since these algebras are abelian). The exact 
sequence relevant for the proof of Levinson's theorem is a non-commutative analog of this, where 
the two commuting coordinate functions in C QO (R) © C 00 (M) are replaced by the operators A and B 
obeying the canonical commutation relations, and then everything is tensorized with the C*-algebra 
of compact operators. 

Let S be the C*-algebra generated by operators of the form f(A) © K and g(B) © K with 
f,g G Co(M) and K in the set JC of compact operators on L 2 (S, z/). Let £ denote the extension of 3 
obtained by allowing / and g to have nonzero finite limits at ±00. Evaluation at infinity of £ gives 
the algebra A which is the subalgebra of (Coo (A) © C QO (B) © C^A) © C^B)) © /C of operators 
having coinciding limits at each of the four corners. This leads to the following short exact sequence 
(see [EE]): 

O^d^Z^A^O, 

It induces a six-term exact sequence in i^-theory |Bla] . leading to a canonical index map Ind : 
Ki(A) — > Ko(3)- As it turns out, using the unitary operators UT where U is defined in equation 
(120]) and J 7 is the Fourier transform, the elements of 3 are mapped into compact operators on £ 2 (Z d ). 
Therefore K (3) = Z |Bla] . Hence the index is an integer. At this level of generality, the index is 
obtained as follows. Given a unitary element S in A, it is lifted to an element Q G £ such that 
ev(fi) = S. Then, ev(f2f2* — 1) = = ev(fi*fi — 1). This means that, if Q can be chosen to be a 
partial isometry, then both QQ* — 1 and Q*Q — 1 are compact projections. Therefore the index of S 
is defined by Ind(S') = Tr(fifi* — Q*Q) which is indeed an integer. At a more computational level, if 
there is a derivation d acting in A and if there is a faithful <9-invariant trace on A, the index of S is 
computed as Ind(S) = Tr(S- 1 9S)/(2vr) pij . 

In the present situation, combining the S- matrix, the operators R± and 1, one for each side of 
the square at infinity, make up a unitary element S tot = S © R + © 1 © R- of A. In addition, in both 
cases, thanks to equation (j42l and to the definition of R± (see Section [331) . this unitary operator 
can be seen as ev(f2_) where f2_ is one of the wave operators. If it can be proved that f2_ G £, then 
Levinson's theorem follows immediately from the following two remarks: 

(i) = 1 and 1 — is the projection onto the pure point spectrum of H. 

(ii) The algebra A admits a trace Tr/i which consists of the trace on /C followed by the integrals over 
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a and h respectively along each side of the square at infinity. In addition, the derivatives 8 a and <9& 
define a derivation 8a acting on A which leaves Tr^ invariant. These derivatives are also given by 
the commutators with A or B respectively. Therefore, the index of S tot is given by 

Ind(5 tot ) = Tr(mr-0*Q) = ^Tta^S^-'OaS^) , 

which is exactly equation (147j) when written in the language used previously. 

The main remaining point is to prove that the wave operators are elements of the C*-algebra £. 
In some sense, most of the technical preparations above have been dedicated to proving precisely 
this point. The formula in Theorem [2] shows that Q± is given by a concatenation of three operators: 
first the operators 0± = J db 0±^ which by Theorem [2] has fibers 0±,& depending continuously on b 
with limits at b = ±00 (due to Proposition [T3|) . then a smooth function of A which also has limits at 
infinity, and finally the imaginary part of the Green matrix also having these properties. One notes 
that the limits in the four corners coincide due to Propositions H2] and [131 For indeed 



lim S b = l — 2\ip±)(ip±\, lim l b = 1 , lim R± a = 1 - 2 \i/j±) (i/j±\ , lim R± a = 1. 

b— >±oo 6— »±oo a— >— 00 a— >+oo 

In conclusion, Q_ G £. □ 



3.9 Example of a point interaction 

Here we discuss the example of the perturbation V = A |0)(0|, A G K, localized on one site. Hence 
n = |0)(0|, A = {0} and, one has L = |A| = 1 even if £ is a polynomial. Furthermore Gq(z) = G (z) 
is a number (and not a matrix of larger size). The behavior of the real and imaginary part of 
Gq(E — zO) can directly be read off Proposition HI Note that, in particular, the imaginary part does 
not vanish on (E_,E + ). Let us introduce the critical coupling constants A± = 1/Gq(E±). Note that 
A_ < and A + > 0. Because the perturbation determinant is 1 — XG (z), the operator H = H + V 
has an eigenvalue smaller than E_ if and only if A < A_, and an eigenvalue larger than E + if and 
only if A > A + . If A = A± there is a threshold singularity at E±. In dimension d = 3 and d = 4 this 
singularity is a threshold resonance, whereas for d > 5 it is a threshold eigenvalue. There is never an 
embedded eigenvalue (also not for polynomial £). 

The scattering matrix given by (|4*4|) differs from the identity only on a one-dimensional subspace: 

9 - (1 TFH 1 I ^-GojE-fQ) 

Hence 

Tr(T E ) = - 8 E m 



% - G (E + zO) 
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From this one readily deduces the winding number Tr(T) if A 7^ A±. For the exceptional cases 
of threshold singularities, we need more precise information about the Green function as given in 
Proposition HI For d = 3 and d > 5 one has 

G (E - i0) = G (E±) + N±(E± - E)+iD±\E - E ± \^ + o{E - E±) . 

The constants satisfy D± > 0, and, for d > 5, one also has N± < 0. Hence, if ip(E) denotes the phase 
of A± x -Gq(E-iO), one has ip(E + ) = and <p{E-) = it with one loop in the clockwise orientation 
for d > 5, while in d = 3 one has <p(E + ) — (p(EJ) = | for A = A±. Hence the total scattering phase is 



Tr(T) = / dE Tr(T E ) = 2n 



for A G (A_, A+) , 

— 1 for A < A_ or A > A + , 

— 1 for A = A± and d > 5 , 

— \ for A = A± and d = 3 . 



This fits with Levinson's theorem. In particular, in the threshold case for d > 5 there is a threshold 
eigenvalue, while for d = 3 there is a threshold resonance (no bound state) but the correction on the 
r.h.s. of f T4Tl) is — \ because dim(«S£ ± ) = 1. 

It is also instructive to complete the argument principle proof in the situation A G [A_, A+]. We 
start from (HSjl and remove 7 as above. Replacing the resolvent identity ( |32|) . then gives 







dz (O|02-# o r 2 |O) f dG 1 



2m A- 1 - G (z) J Go(r) 2m A" 1 - G 



where in the second equality we used (0\(z — H )~ 2 \0) = —d z G (z). Now let us analyze the path 
G (r) using the results of Proposition HI It is always a closed curve going through the real axis twice 
exactly at Gq{EJ) and Gq(E + ). It is positively oriented and the limit curve is approached from the 
inside (as e J, 0). The dimension d now leads to the following crucial differences as to how the real 
axis is crossed. For d = 3 it crosses transversally, as for an circle, while for d > 4, it is a spike 
pointing inward. Hence, if A = A±, the singularity leads to a contribution equal to \ for d = 3 and 
equal to 1 for d > 4. In particular, this shows that Tr(T) = —1 for A = A± also for d = 4, a case 
that was not covered by Theorem [H 

A Boundary values of the Borel transform 

For the convenience of the reader, this section is a reminder of properties of the Borel transform of a 
function defined on the real line. The first result is often called the Plemelj-Privalov theorem |Mus] . 

Lemma 10 Let p : K. — > C be a Holder continuous function of exponent a G (0, 1] with compact 
support. Then, for any (3 such that < (3 < a, its Borel transform 

C, W = f ^ (49) 
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is holomorphic in C \ supp(p) and its boundary value on the real axis is Holder continuous with 
exponent (3. If p is real-valued, then 



G p (E±iO) = Tiwp{E) + f 

Jr 

where f denotes the Cauchy principal value. 



p(e) de 
E-e 



Proof: The holomorphy of G p outside supp(p) is a standard result that will not be proved here. By 
decomposing p into its real and imaginary part, if necessary, there is no loss of generality in assuming 
that p is real- valued. For e > and E eR, G P (E + ie) is given by 

G p (E±ze) = J/e ^E^y + f = R p( E ^) =F *WE,e), 

where 

The first term admits irp(E) as a limit as e | 0. This is because, using the change of variables 
e = E + ex and the Lebesgue dominated convergence theorem, gives 

5? ^' e) = 5? L^Ti p{E + ex) = I^TT = " P{E) - (51) 

Similarly, using the change of variable u = e — E and the symmetry u (-> — u, one obtains 

r°° ?/ dn 

R p (E,e) = J^- 2 ( p (E-u)-p(E + u)). (52) 
Jo u -\- e 

For all a > 0, the part of the integral corresponding to < a < u is also Holder continuous of 
exponent a w.r.t. E, thanks to Lebesgue's dominated convergence theorem. In particular, if E is 
not in the support of p, the integral over u never reaches u = so that R P (E) is Holder continuous 
outside the support of p. However, if p is real valued, R p is the restriction to the complement of the 
support of p (in the real line) of the real part of an holomorphic function and is therefore analytic. 
On the other hand, since p is Holder continuous of exponent a and with compact support, it follows 
that there is a constant K > for which \p(E + 5 ± u) — p(E ± u)\ < K5 a uniformly w.r.t. E and 
u. In particular, \p(E — u) — p(E + u)\ < K(2u) a and 



\p(E + 5 + u) -p(E + u) -p(E + 5-u) +p(E-u)\ < 2Kmin{5 a , (2u) a } < 2 1+a ~^ K5 l3 u a ~^ , 

for any < f3 < a. Using this estimate inside the part of the integral for which u G [0, 1] and thanks 
to the dominated convergence theorem, it follows that lim^o Rp(E, e) exists and is Holder continuous 
of exponent f3 for E G supp(p). The last formula also follows from the above. □ 
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Corollary 2 Let p : R — > C fee k-times differentiable with k-th derivative Holder continuous of 
exponent a G (0, 1] compact support. Then, for any (3 such that < ft < a, its Borel transform 
G p (z) defined by f|4*9l) has a k-times differentiable boundary value on supp(p) for which the kth 
derivative is Holder continuous with exponent (3. 

Proof: Equations ( l5Ti) and ( 1521) show that <9£/ p (e,e) = I p >(E,e) and similarly for R p . Using 
Lemma [TD] gives the result. □ 

Lemma 11 Let p be a real valued C 1 -function on the real line with compact support such that 
p(0) 0. Let denote the partial Borel transform defined by 




Then one has the following. 

(i) G p is holomorphic outside of [0, oo) fl supp(p). 

(ii) The boundary value I p {E) = QmG^(E =p iO) of its imaginary part is ±7cp(E) whenever E > 0, 
it is ±|p(0) at E = and it vanishes for E < 0. 

(iii) The boundary values R P (E) = $t,eG p (E ± zO) is continuously differentiable for E ^ and 
satisfies 

R P {E) = p(0)ln(|£|) + 0(1), asE^O. 

Proof: The proof will use the results given in the proof of Lemma dOj The first claim (i) can be 
proved in the same way indeed. 

(ii) Thanks to equation ( 150|) and with the change of variable e = E + ex, the result of equation 
(15 ip is still valid as long as E > 0. This is because the interval of integration for the variable x is 
[—E/e, oo) in the present case. Then, as e 1 0, the dominated convergence theorem gives the same 
result since this interval of integration eventually becomes R. The same argument applied to E < 
gives and interval of integration [|-E|/e, +oo) which shrink to the empty set as e \. 0. At last, if E = 
the interval of integration is [0, oo), and the result is | instead of n. 

(iii) The argument given in the proof of Lemma fTOl applies also for E ^ 0, so that R p is contin- 
uously differentiable as well. Near E — 0, however, the integral defining R p diverges. This can be 
seen as follows, using equation (joTj) (note the change of sign) 

R P (E,e) = J"de ^ e) _%~ + e l = ~\ [ P(e) d (M(E - e)* + e>)) . 
Since p has compact support, integrating by parts yields 

R p (E, e) = HE" + e 2 )-\ J™ de p\e) hx({E - ef + e 2 ) . 

Note that the second integral converges as e 1 to 
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- / dep'(e) \n(\E-e\) . 
Jo 

which is continuous with respect to E near E = 0. This can be seen by changing e into e — E and 
using the fact that the derivative p' is continuous as well as the dominated convergence theorem. □ 



B Technicalities on certain operator inverses 

Lemma 12 Let A = A* and B > be bounded operators on a Hilbert space. Then Ker(A + iB) = 
Ker(A) D Kei(B) = Ker(A — iB). In addition, the restriction of the operator K = (A + iB)~ l B to 
(Ker(v4) n Kei(B)) ± is well-defined. In particular, ifv G (Ker(v4) HKer(.B))\ thenw = (A + iB)~ l v 
exists as the unique solution w G (Ker(A) D Ker(B)) ± of (A + iB)w = v. 

Proof: If w G Ker(v4) D Ker(B), then clearly (A + iB)w = so that w G Ker(A + iB). Conversely, 
if w G Ker(A + iB), then = + %B)w) = (w\Aw) + i(w\Bw). Since A and B are selfadjoint, 

it follows that (w\Aw) = (w\Bw) = 0. Since B is positive, ||£>u>|| 2 = (w\B 2 w) < \\B\\{w\Bw) = 0. 
Hence w G Ker(B). This also implies Aw = so that w G Ker(A) R Ker(5). The same argument 
shows that Ker(A-tB) = Ker(A)r\Ker(B). Let P be the orthogonal projection onto Ker(74)nKer(5). 
Then it commutes with both A and B. Therefore, both (A + iB) and B restrict to Ran(l — P). 
Moreover, Ran(v4 + iB) = Ker(A — iB) L = Ran(l — P). Therefore the equation (A + iB)w = v 
admits a solution, which is unique modulo Ker(A + iB) which is unique if, in addition, Pw = 0, 
namely if w G (Ker(A) n Ker(i?))" L . □ 

Lemma 13 For an interval I C R, let b G / H> A b G Mat(JV x N, C) and b G / m- B b G Mat(A^ x 
N,C) be real analytic functions satisfying A b = Al and B b > 0. Let P b denote the projection 
on Ran(Sb). Then dim(P&) is constant except for b in a discrete set V C I. We suppose that 
Ker(Afe) n Ker(I?5) is non-trivial only for a discrete set of b 's not lying in V and that those zeros 
are regular singular points for A b , that is, these zeros of b G / h-> A b are of first order. Then 
b G / i—)- K b = (A b + iB b )~ l (B b )?- is well-defined, continuous and real analytic except in those points 
ofV which are not zeros of A b . 

Proof: First of all, by analytic perturbation theory the eigenvalues of B b vary analytically with b 
and therefore dim(Pb) indeed only differs from the almost sure value L on a discrete set V. On V, 
the dimension of P b is then smaller. By the proof of Lemma [T21 if Ker(^4b) D Ker(B b ) is trivial, 
the inverse of A b + iB b exists and, moreover, it is real analytic in b. Thus we may restrict our 
attention to one point, say b = 0, at which Ker(A ) D Ker(B ) is non-vanishing. Then K b is real 
analytic in b G J\{0}. By Riemann's theorem on removable singularities it is sufficient to show 
that K b remains bounded in a neighborhood of b = 0, because then it is already analytic at b = 0. 
Now dim(Ker(5b)) = N — L is constant for 67^0 (possibly after having made / smaller). Again 
by analytic perturbation theory there exists a real analytic function b G / i— > U b of unitaries such 
that U b B b U b = (o °) where rj b > is a diagonal matrix of size L. It is clearly sufficient to show 
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that b G I (->■ U b K b U^ = (U b A b U^ + iU b B b U^) 1 U b (B b ) 1 2U b * is real analytic, which is equivalent to 
supposing that B b is diagonal. Next let us introduce the notations 



u b A b u; 



a b (3 b 
lb 



where a b and 7& are self-adjoint square matrices of size (N — L) x (N — L) and L x L respectively. 
Because Ker(A Q ) D Ker(i? ) is non-trivial, it follows that Ker(a ) H Ker^) is non-trivial as well. 
Moreover, is of size L X (N - L) and (J,J) = C/ 6 P 6 A b (l - P&)£/ 6 *. Furthermore 



u b K b u: 



a b ft \ VO 



where the appearing inverse only exists in the sense of Lemma [121 In order to calculate the inverse, 
let us introduce the Schur complement by 

a b = a b - f3 h {lb + iri b )~ l Pl 

= a b -/3 fe (7 b % 1 7b + ^)~ 1 76% 1 /3b + i PbdbV^lb + Vb)~ l Pt (53) 

which is an (N — L) x (N — L) matrix. The second formula shows also that ^sm(a b ) is non-negative. 
For b 7^ 0, one has Ker(a&) fl Ker(/3^) = {0} so that a b is invertible because then the imaginary part 
is positive. Now the Schur complement formula gives 



Pb \ 1 _ ( o b X ~(y b l Pb{lb + irib] 



-i 



S Ib + iVb) " \-{lb + i?lb) ^jfo 1 (lb + iVb) 1 + {lb + iVb) 1 Pto b 1 (3 b {l b + ir]b) 1 
so that 

uku* (° -Vh^bdb + tVbyHvb)^ 

b b b " Vo (ib + tVbyHvb^ + db + tVby^Zcx^Pbdb + tVb)- 1 ^ 

This shows that U b K b U^ remains bounded in a neighborhood of b = if <J^ 1 f3 b remains bounded, or 
equivalently that F b = cr^ 1 /3 b ('j b r] b ~ 1 'y b + i] b )~^ remains bounded. Setting 

C b = a6-/3fe(7b% 1 7fe + '?6)" 1 7fe% 1 /3b , D b = fi b {^ b % + Vb)~* , 
this can be written as 

F b = {C h + iD h Dt)- 1 D i , 

due to the formula ( 1531) . Now the zero of A b at b = is of first order by assumption, therefore the 
zeros of a b and (3% are of first order as well. It follows that the zero of C b is of first order as well. 
Hence the following lemma is applicable to deduce that F b is bouned. □ 

Lemma 14 For an interval I CR, let b G / h-» C b G Mat(iV x N, C) and b G / M- D b G Mat(iV x 
L,C) be real analytic functions satisfying C b = C£. We suppose that both are invertible except on a 
discrete set of points and that, in case Ker(C&) fl Ker(Z^) is non-trivial for some such b, the zero of 
C b is of first order. Then b G 1 1— > F b = {C b + iD b Dl)^ D b is well-defined and real analytic as well. 
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Proof: Again we focus on one point 6 = for which Ker(Co) D Ker(DQ) is one-dimensional. Then, 
also as above, we choose a basis such that C b is diagonal. In this basis, let us introduce the following 
notations for the matrix entries: 

(bn b 0\ fbrj b f3 b \ 

Cb ~ \ 8 b ) ' D » ~ {bj b a b ) ■ 

Here n b and 5 b are diagonal invertible matrices (for all b) and the size of n b is the multiplicity of the 
zero of C . Then 

n n * ( b 2 (v* b Vb + lllb) K(3 b + 7 *a 6 )\ 

Thus we introduce the Schur complement for (C b + iD b Dl)~ l : 

a b = bn b + i b 2 {i]li] b + 7 6 *7 6 ) - b 2 (n* b (5 b + ^* b a b ) (6 b + i {(3 b (3 b + a* b a b )) \f3 b i] b + a* b lb) ■ 

Note that the appearing is well-defined also at b = because 5q is invertible and the other contribution 
is imaginary. We also note that cr b = bn b + 0(b 2 ). Now we obtain 

m _u , n n* W n - ^ 6 + 1 b2 ^Vb + 7 b *7fc) i Kvth + i* b a b ) V 1 / % 67A 

The upper left corner of the inverse is {<J b )~ l = O^ 1 ) and thus singular, but all the other entries 
remain bounded for b — > because (o"&)~ 1 is always multiplied by a factor b. The singularity of the 
upper left corner, however, is multiplied by b from the second factor. In conclusion, (C b +iD b Dl)~ l D b 
remains bounded and therefore the singularity is removable. □ 
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